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SUMMARY 


The  development  of  field  oriented  techniques  for  determining  the  atmospheric  effects  on  visible  image  transmission 
are  described  in  this  report.  Simplified  analytic  methods  are  presented  for  calculating  the  background  sky,  cloud  and  terrain 
radiance  distribution,  and  in  turn  the  spectral  contrast  transmittance,  of  inclined  paths  of  sight  in  the  atmosphere.  The 
directional  path  radiance  is  determined  as  the  sum  of  the  singly  scattered  sunlight  component  and  a  multiply  scattered 
diffuse  component  that  is  calculated  through  application  of  the  delta-Eddington  approximation.  Performance  tests  of  the 
model  radiance  calculation  vere  carried  out  through  comparisons  with  results  from  more  comprehensive  radiative  transfer 
models  and  by  direct  comparisons  with  high-resolution  scanning  radiometer  measurements  gathered  by  instrumented  air¬ 
craft  in  a  broad  range  of  environmental  conditions.  Techniques  for  the  representation  of  the  optical  parameters  required  as 
model  input  data  are  presented  and  discussed,  including  the  specification  of  cloud  optical  depth  as  a  function  of  cloud  type. 

An  added  feature  of  the  contrast  transmittance  modeling  sysiem  deals  with  the  combined  influences  of  atmospheric 
and  target  factors  in  order  to  calculate  visual  detection  range  for  objects  located  at  the  surface  and  viewed  from  aloft.  Ana¬ 
lytic  representations  of  vision  experiment  data  are  added  along  with  target  information  to  explore  the  relative  sensitivity  of 
detection  range  to  selected  changes  in  target  and  atmospheric  parameters.  Results  are  presented  for  several  examples 
assuming  simplified  targets  and  backgrounds. 
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Frontispiece 

The  analytic  techniques  described  in  this  report  are  a  natural  extension  to  a  broad  program  of  research  into  the  optical  properties  ot  the 
atmosphere  carried  out  by  the  Visibility  Laboratory,  University  of  California.  San  Diego,  in  cooperation  with  and  under  the  sponsorship 
of  the  Air  Force  Geophysics  Laboratory  The  composite  photograph  shown  above,  illustrates  some  elements  of  the  field  program  and 
some  of  the  important  components  of  the  radiance  field  as  measured  by  the  instrumented  aircraft  The  field  site  near  Mt  Rainier.  Wash¬ 
ington,  is  one  of  many  locations  in  the  United  States  and  Europe  where  simultaneous  measurements  of  optical  and  meteorological  vari¬ 
ables  were  gathered  Extensive  high  resolution  measurements  of  the  upper  and  lower  hemisphere  spectral  radiance  fields  for  scenes  such 
as  those  illustrated  by  the  fisheye  camera  photographs  shown  below  (flight  C-359)  were  of  prime  importance  both  for  algorithm  develop¬ 
ment  and  for  validation  of  the  modeling  techniques  over  a  wide  range  of  environmental  conditions. 


ANALYTIC  TECHNIQUES  FOR  ESTIMATING 
IMAGE  TRANSMISSION  PROPERTIES  OF  THE  ATMOSPHERE 

Wayne  S.  Hering 


1.0  INTRODUCTION 

As  set  forth  by  Duntley  (1948),  the  degradation  of 
visible  image  contrast  along  a  viewing  path  of  increasing 
length  is  governed  by  (1)  the  attenuation  of  the  inherent 
background  and  object  radiances  by  air  molecule  and  aero¬ 
sol  particle  scattering  and  absorption,  and  (2)  the  genera¬ 
tion  of  path  radiance  by  molecular  and  aerosol  scattering 
of  incident  light  into  the  path  of  sight  Direct  field  meas¬ 
urements  of  the  basic  optical  properties  required  for  the 
determination  of  these  factors  for  selected  path  segments 
cannot  be  made  simply  and  inexpensively.  Thus,  opera¬ 
tional  estimates  of  contrast  transmittance  must  be  derived 
from  mode!  simulations  that  relate  the  aerosol  directional 
scattering  and  absorption  properties  to  generally  available 
meteorological  observations  and  forecasts  (e.g.  Shettle  and 
Fenn,  1979)  Additional  problems  stem  from  the  com¬ 
plexity  of  multiple  scattering  processes  and  the  relatively 
large  amount  of  computer  time  and  capacity  required  for 
treating  the  radiative  transfer  calculations  completely  as  a 
physical  problem  Effective  simplifications  must  be  made 
to  yield  fast  yet  consistent  calculations  of  the  radiance  dis¬ 
tribution.  The  successful  development  of  techniques  for 
real  time  estimates  of  the  atmospheric  effects  on  visible 
image  transmission  requires  a  reasonable  balance  between 
the  desired  accuracy,  the  complexity  of  the  modeling  pro¬ 
cedures,  and  the  reliability  and  representativeness  of  the 
available  input  data.  Since  these  factors  may  vary  consid¬ 
erably  from  one  application  to  the  next,  it  is  important  to 
retain  reasonably  good  flexibility  and  completeness  in  the 
analytic  approximations  used  for  the  estimates. 

The  development  and  testing  of  field  oriented  tech¬ 
niques  for  calculating  directional  path  radiance  and  the 
object-background  contrast  transmittance  along  inclined 
paths  of  sight  in  clear-hazy  atmospheres  were  described  in 
a  previous  report  (Hering,  1981).  This  paper  addresses 
extensions  and  refinements  to  the  modeling  procedures 
and  their  validation.  In  particular,  techniques  are 
presented  for  (1)  estimating  the  background  radiance  and 
visible  image  transmission  characteristics  in  overcast  sky 
conditions  from  conventional  observations  or  forecasts  of 
cloud  type  and  altitude.  As  in  the  case  of  prior  studies, 
technique  development  was  based  extensively  on  experi¬ 
mental  data,  including  a  broad  series  of  simultaneous  opti¬ 
cal  and  meteorological  measurements  gathered  with  a  spe¬ 
cially  instrumented  aircraft  The  high  resolution  profile 
measurements  were  obtained  over  a  period  of  several 
years  in  a  wide  range  of  meteorological  and  geographical 


conditions  in  the  United  States  and  western  Europe. 

Further  extensions  of  the  modeling  concepts  deal 
with  the  combined  effects  of  both  environmental  and  tar¬ 
get  factors  on  visual  detection  range  Although  indepen¬ 
dent  assessment  can  be  made  of  the  changes  in  contrast 
transmittance  that  are  associated  with  natural  changes  in 
the  relevant  atmospheric  variables,  the  impact  of  these 
changes  on  target  detection  distance  depends  markedly  on 
such  factors  as  inherent  target  contrast,  target  size,  and 
glimpse  time  Analytic  representations  of  visi-  n  experi¬ 
ment  data  (Blackwell,  1946  and  Taylor.  I9M1  were  added 
to  the  contrast  transmittance  model  to  explore  the  relative 
sensitivity  of  visual  detection  distance  to  ■  selected  range 
of  target  and  atmospheric  variables  Results  are  presented 
for  several  examples  assuming  simplified  largets  and  back¬ 
grounds 

A  review  of  the  basic  equations  and  a  summary  of 
the  analytic  approximations  used  for  calculating  path  radi¬ 
ance  and  for  estimating  visible  spectrum  contrast  transmit¬ 
tance  in  cloudless  atmospheres  are  given  in  Section  2. 
Also  presented  in  Section  2.  are  comparative  radiative 
transfer  calculations  between  the  simplified  model  and 
more  complete  computational  methods.  In  Section  3, 
techniques  are  presented  for  the  representation  of  the 
input  variables  that  are  required  for  the  calculation  of 
background  radiance  and  spectral  contrast  transmittance. 
Extension  of  the  modeling  concepts  to  deal  with  image 
transmission  characteristics  for  viewing  paths  below  an 
overcast  cloud  layer  are  described  in  Section  4.  Results  of 
validation  tests  comparing  model  calculations  with  experi¬ 
mental  measurements  of  the  directional  radiance  distribu¬ 
tion  gathered  by  instrumented  aircraft  are  given  in  Sec¬ 
tion  5.  The  results  of  sensitivity  analyses  which  give  a 
measure  of  the  relative  importance  of  individual  target  and 
environmental  factors  on  the  determination  of  target 
detection  range  are  summarized  in  Section  6. 

2.0  THEORY  AND  ANALYTIC  APPROXIMATIONS 

Although  some  refinements  and  additions  to 
doudless-sky  modeling  procedures  are  presented,  this  sec¬ 
tion  is  primarily  a  review  of  theory  and  analytic  approxi¬ 
mations  given  in  Hering  (1981) 

Neglecting  turbulence  effects,  the  equation  for  the 
apparent  spectral  radiance  of  the  background  h  at  range  r 
along  the  path  of  sight  specified  by  zenith  angle  6  and 
azimuthal  angle  d>,  can  be  written  (Duntley  et  al.  1957) 
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„  4,  (;,»,<*)-  kLo(z,,0.  0.  <t>>  .  0) 

where  »£.„  is  the  inherent  background  radiance  at  target 
altitude  z, , 

T,  -  exp  -/  a  (r ) dr 


is  the  path  transmittance,  a  is  the  volume  attenuation 
coefficient,  and  L ’  is  the  path  radiance  produced  by  the 
scattering  of  light  from  the  sun  and  from  the  surrounding 
sky,  clouds  and  terrain  into  the  path  of  sight. 

The  path  radiance  is  given  by 

T 

L*(Z'9,  <t>)  -  f  L.  6.  <b)T,  (z',0)dr'  .  (2) 


where  L.U',  0,6)  is  the  path  function,  defined  as  the 
point  function  component  of  path  radiance  generated  by 
the  directional  scattering  of  light  reaching  that  point  of  the 
path.  The  expression  for  the  path  function  can  be  written 
in  terms  of  its  contributions  from  the  scattering  of  solar 
(lunar)  scalar  irradiance  ,<  (r )  and  from  the  scattering  of 
sky  and  earth  radiances  L(z,0\ 6')  as  follows  (eg.  Gor¬ 
don,  1969b): 

4,(z,0,  <b)  -  ,((z)<r(z )  +  f  L  (z  fi'.<b')(r(z  .  (3) 


where  a(z  ,p)  is  the  directional  volume  scattering  function 
at  angle  /3  between  the  path  of  sight  and  direction  of  the 
source  light,  dfl  is  an  element  of  solid  angle,  and  /3S  is 
the  scattering  angle  with  respect  to  the  sun. 

In  turn,  the  spectral  contrast  transmittance  of  the 
path  of  sight  can  be  expressed  (see  Duntley  et  at.  1957) 
directly  as  the  product  of  the  path  transmittance,  Tr(z,8), 
and  the  ratio  of  the  inherent,  bL0,  and  apparent,  bL, , 
background  radiances  as  follows: 

C,(r .#.4)/ C„(z,,@.4>)  -  T,(z,e)l,L<)(zl,e,<t>)/ bL,(z,8,<t>M) 

where  C,  -  (,L,  -  »/*)/ bL,  is  the  apparent  target  con¬ 
trast  at  path  length  r,  C« -(,!»- *£»)/*£»  is  the 
inherent  target  contrast  at  altitude  z,  and  -Lj  is  the 
inherent  target  radiance.  Thus,  the  contrast  transmittance 
of  a  given  path  does  not  depend  upon  intrinsic  target 
characteristics  but  is  a  function  only  of  the  directional 
radiance  distribution  in  the  atmosphere  and  the  path 
transmittance.  The  expression  is  strictly  applicable  only 
for  monochromatic  radiation  but  may  be  applied  with  good 
approximation  to  reasonably  broad  spectral  bands  in  the 
visible  portion  of  the  spectrum. 

From  Eq.  (3)  we  note  that  the  path  function  and  in 
tum  the  contrast  transmittance  as  calculated  from  Eq  (4) 
depend  upon  the  direction  of  the  viewing  path  relative  to 


the  distribution  of  light  reaching  the  path.  The  first  term 
on  the  right  hand  side  of  Eq.  (3)  is  the  contribution  of  pri¬ 
mary  scattering  of  direct  solar  irradiance  The  directional 
volume  scattering  function  may  be  expressed 

ir(z.ii)  =  P(:.fi>s<zl  .  (5) 

where  s(z)  is  the  total  volume  scattering  coefficient,  and 
P(z,fS )  is  the  single  scattering  phase  function  for  com¬ 
bined  Rayleigh  plus  aerosol  particle  scattering  (see  discus¬ 
sion  in  Section  3.3).  which  defines  the  probability  that 
incident  radiation  will  be  scattered  in  the  direction  given 
by  scattering  angle  /3.  The  phase  function  varies 
significantly  with  the  scattering  properties  of  the  aerosol 
particle  distribution  in  the  atmosphere 

Substituting  Eq  (5)  in  Eq  (3),  the  expression  for 
L.izjKd )  becomes 

)  -  ,t  (r  )P{z  ,p)s  (:  l 

L(:.H'.6')P(:.0  .dll'  (6) 

Modelling  techniques  for  the  calcula'  of  the  diffuse 

component  of  the  path  function  are  c  *  in  the  fol¬ 

lowing  paragraphs. 

2.1  Calculation  of  the  Path  Function 

The  second  term  on  the  right  hand  side  of  fcq  <6>  is 
the  component  of  the  path  function  resulting  from  the 
scattering  of  diffuse  radiance  reaching  the  path  from  the 
surrounding  sky  and  terrain  It  has  a  directional  depen¬ 
dence  due  to  the  asymmetry  in  the  background  sky  and 
earth  radiance  distribution  and  the  scat,  mg  phase  func¬ 
tion  asymmetry.  Precise  numerical  calculation  of  the  nath 
radiance  resulting  from  the  complex  multiple  scattering 
processes  requires  large  amounts  of  computer  time  For 
this  reason,  rapid  appro-  imate  methods  are  employed 
extensively  for  radiation  transfer  calculations  The 
appropriate  choice  of  computational  method  from  among 
the  variety  of  available  methods  depends  upon  the  results 
desired  for  the  application  at  hand. 

While  it  is  important  to  retain  complete  directional¬ 
ity  for  calculation  of  the  path  radiance  component  due  to 
single  scattering  of  direct  solar  radiance,  approximate 
hemispherical  two  stream  methods  can  be  used  effectively 
for  fast  calculation  of  the  path  radiance  component  due  to 
scattering  of  the  background  sky  and  terrain  radiances  at 
any  point  and  direction  provided  that  the  asymmetric 
influence  of  the  prominent  forward  scatter  peak  and  the 
irradiance  profile  are  managed  adequately.  The  delta- 
Eddington  approximation  introduced  by  Joseph.  Wiscombe 
and  Weinman  (1976)  satisfies  the  requirement.  It  extends 
the  standard  Eddington  approximation,  which  assumes  a 
simple  cosine  dependence  of  the  single  scattering  phase 
function,  through  approximation  of  the  phase  function, 
Pj,  by  a  truncated  forward  scatter  peak  and  a  two-term 
phase  function  expansion. 
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4ir  Pd(0)  -  2/‘Ml  -  cos/9)  +  (I  -  /'Ml  +  3g'cos0)  (7) 

where  /'  is  the  fractional  scattering  represented  by  the  for¬ 
ward  peak  and  g'  is  the  asymmetry  factor  of  the  truncated 
phase  function  In  effect,  the  delta-Eddington  approxima¬ 
tion  transforms  most  of  the  enhanced  radiance  in  the  solar 
aureole  into  the  direct  solar  flux  component,  and  it  is  gen¬ 
erally  assumed  that 

f'{:  1  -  g:{z)  .  (8) 

This  assumption  is  commensurate  with  representation  of 
the  actual  phase  function  with  a  single  term  Henyey- 
Greenstein  phase  function  Alternate  expressions  are 
required  (see  McKellar  and  Box.  1981)  for  other  phase 
function  representations  such  as  the  2-term  Henyey- 
Greenstein  functions.  Joseph  et  at.  (1976)  show  that  cal¬ 
culations  of  radiative  transfer  with  the  delta-Eddington 
approximation  can  be  carried  out  with  the  standard 
Eddington  computer  code  (Shettle  and  Weinman,  1970) 
with  the  following  changes  of  variable 

At  '  -  ( 1  -  <uf‘)  At  ,  (9) 


<u’  -  Old  -  /’)  /  (l  -  wf )  .  (10) 

and 

<l-/'>  .  (U) 

where  At  is  the  optical  thickness  of  the  layer  and  <a  is  the 
single  scattering  albedo. 

As  an  integral  part  of  the  technique  for  estimating 
directional  contrast  transmittance,  the  products  of  the 
Eddington  computer  program  (Shettle  and  Weinman, 
1970)  are  used  directly  to  calculate  the  second  term  on  the 
right  hand  side  of  Eq.  (6). 

From  the  standard  Eddington  approximation,  the 
diffuse  radiance  is  assumed  to  be  given  by 

L  (z  -  Ld  Or )  +  L0  (z  Icosfl'  .  (12) 

As  shown  by  Shettle  (1981),  if  we  substitute  Eq.  (12)  and 
the  delta-Eddington  approximations  given  by  Eq.  (7),  (8) 
and  (11)  into  the  last  term  of  Eq  16)  and  integrate  overfl' 
and  <6',  we  have 

L(z,0'.<b,)P(z.(i)s  (z)dn 


-  s  (z )  \Ld  (z )  +  gL0  (z  )cos@!  .  (13) 

A  recent  refinement  in  the  above  expression  for  the 
diffuse  component  of  the  path  function  consists  of  an 


additional  term  to  help  account  for  the  effects  of  the 
azimuthal  asymmetry  in  the  background  sky-terrain  radi¬ 
ance  distribution.  The  term  is  introduced  as  a  rough 
approximation  to  an  iteration  of  the  singly  scattered  sun¬ 
light  component  of  the  path  function.  Its  magnitude  is 
assumed  proportional  to  the  single  scattering  component, 
,tiz)P(z,p),  and  proportional  to  the  fractional  contribution 
of  the  total  diffuse  energy,  4n  LD(z ),  to  the  sum  of  the 
total  diffuse  energy  and  singly  scattered  sunlight  com¬ 
ponent  at  the  level  of  computation.  In  other  respects,  the 
modification  conforms  with  the  extension  of  the  delta- 
Eddington  model  introduced  by  Davies  (1980).  in  particu¬ 
lar,  the  contribution  of  the  asymmetry  term  is  proportional 
to  sme.sinft,  and  cos*.  The  modified  expression  for  the 
diffuse  component  of  the  path  function  is 

L  Iz.tf'.  <b'lP(z  .«•• '.vii 


-  s  (2 )  I Ld  (z )  +  gLj)  (z  IcosS  +  gLp  (:  > cos d> )  (14) 

where 

4rr  iD(z)Ptz,0)  sf  (z  Isinft,  smi- 

L°Az ' -  MrrDTrnrpTz7f) tt- >"  ■  n 5) 

Commensurate  with  the  delta-Eddington  approximations 
and  Eq.  (9),  the  expression  for  the  solar  scalar  irradiance 
at  altitude  z  is 

s(  (z )  *  jt  (oolexpt-T,  /  cos#, )  .  (16) 

where  ,«(<»)  is  the  extraterrestrial  solar  scalar  irradiance. 
The  optical  depth  rx  is  assumed  equal  to  the  delta- 
Eddington  optical  depth,  r',  everywhere  except  within  the 
forward  peak  (/?  <  25  deg)  where  rx  is  equal  to  the 
unmodified  value  of  r  Substituting  Eq s.  (13)  and  114) 
into  Eq.  (6),  the  expression  for  the  path  function  becomes 

L,(z  ,0,  A1  -  s  (z)  [P(z  ,0 j  ),<  (°°  )exp(-  t„  i  cost/  ) 

+  Lp(z)  +  gLDtz)eosf)  +  gLp  (z)cos4>l  (17) 

2.2  Calculation  of  Sky  and  Terrain  Radiance 

For  an  assumed  plane  parallel  and  horizontally 
homogeneous  atmosphere,  the  directional  path  radiance, 
L’{z  ,0.  ®),  and  the  radiance,  L  (z  ,6,  <t>),  can  now  be  cal¬ 
culated  from  Eqs.  (1),  (2)  and  (17)  through  finite  summa¬ 
tion  over  adjacent  atmospheric  layers  using  the  trapezoidal 
rule.  For  upward  paths  of  sight,  the  inherent  background 
radiance  of  clear  sky  at  target  altitude  is  given  by 

^(zj.e,*)-  Ll  (z,.«.<*)  .  (18) 

where  LL  (z,  ,9,<b)  is  the  path  radiance  as  determined  at 


altitude  z,  for  the  slant  path  from  the  top  of  the  atmo¬ 
sphere  to  z, .  Looking  downward,  the  inherent  background 
radiance  for  the  case  of  uniform  Lambertian  reflectance 
and  a  horizontal  surface  is  given  by 

4Z.„(2,,e,<*)  -  Tn(z,, 8)R{9.  <t>)E(0,d)hr 

+  L’b,,e,<b)  ,  (19) 

where  RW,<t>)  is  the  local  surface  reflectance, 
£(0,<f)  —  *,  (O)cos0,  +  rr[£p(0)  +  2/3L0 ■  (0) ]  is  the 
downwelling  irradiance  at  the  surface,  and  Tn  is  the 
transmittance  of  the  slant  path  from  the  surface  to  target 
altitude  z, . 

2.3  Comparisons  with  Other  Model  Calculations 

Tests  were  carried  out  to  evaluate  the  performance 
of  the  simplified  model  for  the  radiance  calculations  as 
described  above.  Model  calculations  of  radiance  distribu¬ 
tions  were  compared  with  the  results  calculated  from  more 
comprehensive  radiative  transfer  models,  which  are 
mathematically  precise  but  in  general  require  relatively 
large  computer  capacity  and  lime  for  the  radiative  transfer 
calculations.  For  these  tests,  the  comparative  accuracy 


was  determined  under  the  assumption  that  the  physical 
system  is  perfectly  observed,  and  the  input  data  are  identi¬ 
cal  for  the  approximate  operational  model  and  the 
comprehensive  numerical  models.  The  set  of  examples 
used  for  the  analysis  were  selected  from  those  proposed  by 
the  ad  hoc  Working  Group  of  the  Radiation  Commission. 
International  Association  of  Meteorology  and  Atmospheric 

Table  1.  Summary  of  numerical  comparison  cases  lisied 
in  the  Working  Group  Report  (Lenoble.  1977)  of  the 
Radiation  Commission,  IA.MAP,  for  5  reference  atmo¬ 
spheres 
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fig.  1.  Comparative  radiance  calculations  for  optical  depth  0.1  for  reference  cases  1,2  and  3  Model  results  (dashed  lines)  are  compared 
with  the  spherical  harmonics  method,  (Lenoble,  1977),  (solid  lines).  Optical  parameter  data  for  each  case  are  given  in  Table  1 


Physics  (Lenoble,  1977,  see  Table  1).  The  results  of  the 
spherical  harmonic  method  listed  in  the  Working  Group 
report  were  chosen  for  a  reference  standard  since  com¬ 
parative  calculations  using  this  method  were  presented  for 
all  five  test  cases.  A  single  layer,  plane  parallel  and  homo¬ 
geneous  atmosphere  with  a  black  underlying  surface  was 
prescribed  in  each  case.  Three  cases  have  a  Haze  L  (Deir- 
mendjian,  1969)  aerosol  particle  distribution  with  an  opti¬ 
cal  thickness  of  1,  and  the  other  two  have  assumed  cloud 
particle  distributions  with  an  optical  thickness  of  64.  The 
phase  function  for  single  scattering  and  the  single  scatter¬ 
ing  albedo  are  specified  for  each  case  Comparative  results 
of  the  radiance  calculations  are  illustrated  in  Figs.  1-4. 

It  is  important  to  note  that  the  agreement  among 
the  more  rigorous  numerical  methods,  as  given  in  the 
Working  Group  report,  is  generally  within  a  few  percent 
depending  upon  the  completeness  of  scattering  history 
retained  in  the  calculations.  We  see  from  Fig.  1  that  the 
departures  of  the  simplified  model  radiance  calculations 
from  the  reference  values  for  the  haze  eases  are 
significantly  larger,  yet  the  indicated  accuracy  is  good  when 
considered  in  context  with  the  uncertainties  and  complete¬ 
ness  of  the  observational  data  base  from  which  the  esti¬ 
mates  must  be  uerived  in  operational  practice.  The 


simplified  model  requires  only  a  few  seconds  of  large  scale 
computer  time  for  each  case,  and  the  model  is  designed 
for  easy  microprocessor  application  for  real  time  require¬ 
ments.  The  approximations  inherent  in  the  computation¬ 
ally  fast  model  serve  to  smooth  out  the  asymmetry  in  the 
calculated  radiance  component  due  to  the  scattering  of 
incident  diffuse  light;  hence  it  is  expected  that  the  error  of 
estimate  would  continue  to  increase  somewhat  in  the  solar 
azimuthal  plane  {<4>  -  0-180  deg)  as  the  solar  zenith  angle 
increases  beyond  60°. 

The  relative  effectiveness  of  the  azimuthal  asym¬ 
metry  term,  g  Ly cos<6,  is  illustrated  in  Table  2.  Model 
calculations  for  Case  3  (6,-60  deg,  o>— 0.9)  are  listed 
with  and  without  the  azimuthal  component  contribution 
for  comparison  with  the  reference  model  results.  It  is  evi¬ 
dent  that  the  accuracy  of  Eq.  (17)  is  much  better  with  the 
azimuthal  correction.  However,  systematic  errors  remain 
such  that  further  study  of  this  aspect  of  the  modeling  pro¬ 
cedures  is  particularly  important. 

The  relative  performance  of  the  model  calculations 
of  radiative  transfer  in  a  high  extinction  environment  is 
illustrated  in  Figs.  3  and  4.  Both  of  the  cloud  examples  in 
the  reference  series  assume  an  overhead  sun  position. 
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Table  2.  Comparison  of  Case  3  radiance  distributions  as  calculated  by  the  Spherical  Harmonics  method  Icolumn  SH).  the  simplified  model 
without  azimuthal  asymmetry  term  (column  A),  and  simplified  model  with  the  azimuthal  asymmetry  term  Icolumn  B)  The  solar  zenith 
angle  is  60  deg.  Calculations  are  shown  for  a  level  near  the  top  of  the  haze  layer  (t  -  0.1)  and  a  level  near  the  base  (t  -  0  75)  Note  that 
the  azimuthal  term,  g  LD-  cos*.  is  zero  for  <6-90  deg,  hence  columns  A  and  B  are  identical 
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fig.  3.  Same  a*  Fig.  I  except  for  cloud  case  4  at  optical  depths 
3.2,  12.8,  and  48  and  no  absorption. 


Fig.  4.  Same  as  Fig.  3  except  for  cloud  case  5  with  single  scatter¬ 
ing  albedo  of  0.9. 
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Tne  agreement  with  the  reference  model  calculations 
holds  well  for  both  cloud  cases,  although  the  disparity 
increases  with  increasing  depth  for  the  case  witlt  strong 
absorption.  The  simplified  model  calculations  of  in-cloud 
radiance  were  carried  out  through  division  into  sub-layers 
with  boundaries  corresponding  to  levels  where  reference 
model  calculations  were  available.  These  subdivisions 
represent  large  increments  of  optical  depth  in  the  cloud 
cases  However,  the  estimates  if  in-cloud  radiance  are 
good  in  general  over  the  full  range  of  optical  depth  and 
observing  angles  shown  in  Fig.  3  and  4.  A  notable  excep¬ 
tion  is  for  the  data  point  directly  upsun,  0=tfs=d>= 0,  at 
optical  depth  3.2  The  calculated  values  of  in-cloud  radi¬ 
ance  for  this  data  point  result  in  overestimates  by  at  least 
an  order  of  magnitude  and  plot  offscale  in  Figs.  3  and  4. 
The  modeling  procedures  require  important  refinement  in 
layer  by  layer  resolution  in  order  to  handle  combined 
effects  of  the  rapid  attenuation  of  the  direct  solar  beam  in 
the  cloud  near  the  upper  limit  of  the  cloud  layer  and  the 
strong  change  in  the  single  scattering  phase  function  with 
scattering  angle  for  /3  less  than  5  deg.  The  large  errors  are 
confined  to  viewing  paths  within  a  few  degrees  of  the  sun 
and  to  the  region  near  upper  limit  of  the  cloud  layer.  In 
this  region,  the  singly  scattered  sunlight  component  makes 
a  significant  contribution  to  the  path  function  and  it 
changes  rapidly  along  the  path. 

3.0  APPLICATION  OF  ANALYTIC  TECHNIQUES 

Environmental  data  and  forecasts  relevant  to  the 
determination  of  the  background  spectral  radiance  distri¬ 
bution  and  image  transmission  characteristics  may  be 
available  in  many  forms  and  with  varying  degrees  of  com¬ 
pleteness  and  representativeness.  Emphasis  has  been 
placed  on  the  development  of  a  general  modular  format 
for  requisite  data  entry  seeking  to  take  full  advantage  of 
all  available  information  important  for  a  particular  applica¬ 
tion  of  the  radiance  model  A  summary  of  input  variables 
and  alternate  techniques  for  the  specification  of  the  atmos¬ 
pheric  properties  are  included  in  Appendix  A.  Basic  data 
entries  include  the  wavelength  representative  of  the  sensor 
spectral  characteristics,  solar  zenith  angle,  and  the  obser¬ 
vational  paths  of  interest  for  the  problem  at  hand. 
Specification  of  the  underlying  surface  reflectivity  is  dis¬ 
cussed  below  in  Section  3.1. 

Several  options  are  available  for  the  input  of  specific 
atmospheric  variables.  One  may  introduce  as  many 
atmospheric  layers  as  warranted  by  the  accuracy  and  com¬ 
pleteness  of  the  meteorological  observations  or  forecasts 
available  for  the  specification  of  the  optical  parameter 
profiles  Minimum  information  for  each  designated  layer 
includes  the  altitude  limits,  the  average  scattering 
coefficient,  the  average  absorption  coefficient  and  the  sin¬ 
gle  scattering  phase  function  for  aerosol  particle  scattering. 
Techniques  for  the  representation  of  these  variables  are 
discussed  in  Sections  3.2  and  3.3. 

3.1  Representation  of  Surface  Reflectance 

The  background  surface  reflectance  contributes  to 
the  calculated  radiance  distribution  in  two  ways.  First,  the 


general  or  broad  area  average  surface  reflectance  is 
entered  to  calculate  the  component  of  path  radiance  that  is 
generated  by  light  reaching  the  path  through  scattering 
from  the  underlying  surface  and  in  turn  scattered  in  the 
direction  of  the  sensor  as  given  by  Eq.  (6).  For  a  land 
surface,  the  reflectance  is  assumed  to  conform  with 
Lambert’s  law,  so  that  the  resultant  radiance  is  indepen¬ 
dent  of  observation  angle  and  depends  only  on  the 
downwelling  irradiance  and  the  average  surface 
reflectance. 

Second,  the  local  background  reflectance,  which  may 
differ  from  the  area  averaged  reflectance  discussed  above, 
is  entered  for  the  determination  of  the  inherent  back¬ 
ground  radiance  and  its  contribution  to  the  apparent  spec¬ 
tral  radiance  as  given  by  Eqs.  (1)  and  (19). 

The  above  algorithms  deal  with  surface  reflectance 
over  land  areas.  Another  option,  for  use  primarily  over 
an  open  water  surface,  assumes  specular  surface  reflection. 
Here  the  surface  radiance  is  a  function  of  both  the  viewing 
angle  and  the  downwelling  radiance  distribution  For  the 
water  option,  the  Fresnel  reflectance  of  the  water  surface 
as  a  function  of  observation  angle  is  set  equal  to  average 
values  for  a  sea  surface  roughened  by  light  surface  winds 
(4m  see' 1 ) .  These  values  as  calculated  by  Gordon 
(1969a),  assume  that  the  crosswind  and  upwind  wave 
slope  probabilities  of  Cox  and  Munk  (1954)  may  be 
approximated  by  a  single  circular  distribution  of  wave 
slopes. 

3.2  Representation  of  Total  Volume 
Scattering  Coefficient 

The  spatial  distribution  of  total  volume  scattering 
coefficient  deserves  prime  consideration  since  it  is  the 
major  determinant  of  visible  spectrum  contrast  transmit¬ 
tance.  Techniques  for  the  specification  of  the  scattering 
coefficient  profile  have  been  investigated  during  the  course 
of  the  aircraft  measurement  and  analysis  program.  Some 
of  the  results  of  these  studies  and  their  application  to 
operational  modeling  procedures  were  summarized  by 
Hering  (1981).  A  brief  review  is  given  here  for  immedi¬ 
ate  reference. 

For  profile  modeling  purposes,  it  is  important  to 
consider  a  conservative  measure  of  scattering  coefficient 
that  in  the  absence  of  local  aerosol  particle  sources  or 
sinks  does  not  change  appreciably  following  the  air 
motion.  The  optical  scattering  mixing  ratio,  (?(z),  is  such 
a  parameter.  As  the  vertical  mixing  within  an  identifiable 
atmospheric  layer  becomes  more  complete,  Qiz)  becomes 
more  constant  with  height  within  the  layer. 

The  optical  scattering  ratio  is  defined 

Q(z)  -  s  (z )/  s*  (z)  ,  (20) 

where  s*(z)  is  the  total  volume  coefficient  for  Rayleigh 
scattering  at  altitude  z .  Note  that 

s(r)  -  s„  (r)  +  sM(r)  -  o>(z)a(z)  ,  (21) 
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where  a(z)  is  the  extinction  coefficient  and  sM(z)  is  the 
aerosol  scattering  coefficient.  It  follows  that  the  aerosol 
scattering  ratio  is  given  by 

sv(r)/j*(z)-  Q(z)~  1  .  (22) 

The  aerosol  scattering  ratio  also  would  be  constant  under 
conditions  of  complete  aerosol  mixing.  An  additional 
computational  advantage  of  a  scattering  ratio  representa¬ 
tion  is  that  it  provides  normalization  with  respect  to  both 
density  altitude  and  wavelength. 

Profiles  of  Q(z)  derived  from  the  extensive  series 
of  airborne  optical  measurements  made  by  the  Visibility 
Laboratory,  reveal  large  variability  depending  upon  the 
aerosol  particle  source  strength  and  the  nature  of  the  con¬ 
vective  and  turbulent  mixing  processes  The  problem  is  to 
model  the  essential  characteristics  of  the  Q(z)  profiles  in  a 
v,  jy  that  recognizes  operational  observing  and  forecasting 
limitations  yet  takes  maximum  advantage  of  existing  capa¬ 
bilities.  A  prominent  feature  of  the  daytime  aircraft 
soundings  over  inland  areas  was  the  marked  tendency  for 
Q(z)  to  remain  essentially  constant  with  height  within  the 
boundary  layer,  and  also  in  the  relatively  haze-free  region 
of  the  upper  troposphere  above  the  primary  haze  layer  It 
should  be  emphasized  that  the  assumption  of  constant 
scattering  ratio  with  height  does  not  hold  well  for  ground- 
based  stable  layers  with  little  vertical  mixing  such  as  those 
associated  with  the  nocturnal  formation  of  fog  However, 
for  application  to  problems  of  contrast  transmittance  in 
hazy  atmospheres  in  the  daytime  following  the  dispersion 
of  any  surface  inversion  existing  at  sunrise,  a  simple  3- 
layer  troposphere  model  with  constant  Q(z)  in  each  layer 
provides  in  most  cases  a  good  first  approximation  of  the 
aerosol  scattering  profile  Thus,  the  forecasting  problem  is 
reduced  to  the  prediction  of  the  upper  altitude  limit  of 
mixed  boundary  layer  and  the  avenge  scattering  ratio 


within  each  layer.  An  evaluation  of  model  representations 
of  the  high  resolution  scattering  ratio  profiles  as  measured 
during  several  deployments  of  the  instrumented  aircraft  in 
western  Europe  is  given  in  Hering  (1981). 

Three  examples  of  scattering  coefficient  profiles  as 
measured  by  the  integrating  nephelometer  are  shown  in 
Fig.  5.  Each  experimental  flight  usually  included  profile 
measurements  with  4  spectral  filters  having  peak 
wavelengths  of  475,  550,  660,  and  750  nm.  The  4  profiles 
were  measured  sequentially  over  approximately  a  2-hour 
period.  Data  for  individual  profiles  were  extracted  at  30-m 
intervals.  Figure  5  illustrates  some  recurrent  features  of 
the  lower  tropospheric  profiles;  the  well  mixed  boundary 
layer  in  late  morning  and  afternoon  over  inland  areas  of 
northern  Europe  (Soesterberg  profile),  the  multi- 
structured  profile  in  low-level  fog  conditions  offshore  near 
Rodby,  Denmark,  and  the  deep  haze  layer  typical  of  the 
summer  profiles  measured  in  the  Mediterranean  Sea  area 
(Trapani,  Sicily). 

Cursory  analysis  of  a  large  number  of  simultaneous 
aircraft  measurements  of  temperature  and  scattering 
coefficient  profile  structure  indicate  that  a  temperature 
lapse  rate4 equal  to  the  standard  lapse  rate  of  6.5°C  per  km 
or  greater  will  support  rather  complete  vertical  mixing  of 
aerosol  particles  in  the  size  range  contributing  predom¬ 
inantly  to  extinction  in  the  visible  portion  of  the  spectrum 
On  the  other  hand,  atmospheric  layers  of  thickness  greater 
than  about  300  m  and  lapse  rate  less  than  about  4.5°C  per 
km  significantly  inhibit  vertical  mixing  and  in  general 
appear  sufficient  to  cap  the  mixed  boundary  layer. 

3.3  Representation  of  Single  Scattering 
Phase  Function  and  Single  Scattering  Albedo 

The  single  scattering  phase  function  P(zji)  as 
employed  in  Eq.  (6)  is  a  bulk  parameter  of  the  atmos- 


OPTICAL  SCATTERING  RATIO,  Q(Z> 


Fig.  5.  Scattering  ratio  profiles  as  measured  for  (light  C-372  (Soesterberg,  Netherlands),  (light  C-390  (Rodby,  Denmark)  and  (light  C-462 
(Trapani,  Sicily),  The  heavy  dashed  line  is  an  objective  model  representation  of  the  profile.  The  spectral  nephelometer  measurements 
were  corrected  for  the  effects  of  incomplete  purging  of  the  instrument. 


-  8  - 


pheric  layer,  representing  the  combined  aerosol  and  Ray¬ 
leigh  phase  functions.  It  is  given  by 


/*0,z)  - 


Pg(fi)  +  (Q(z)~  l I4M 


where  Ph(&j)  is  the  phase  function  for  aerosol  scattering 
and  the  theoretical  Rayleigh  phase  function  is 


So  that  the  fraction  of  radiation  absorbed  for  each  single 
photon  collision  is  1  -  u>(z).  Commensurate  with  other 
variables,  <u(z)  is  introduced  as  a  bulk  parameter  and  is  a 
measure  of  the  combined  absorption  effects  of  both  air 
molecules  and  aerosol  particles.  In  the  visible  spectrum, 
biCz)  in  the  boundary  layer  ranges  from  1.0  (no  absorp¬ 
tion)  to  less  than  0.7  in  urban  atmospheres  with  appreci¬ 
able  carbon  concentration. 


Pr  (£)  “  3(1  +  cos2/})  /  16rr  .  (24) 

The  phase  function  for  single  scattering  has  the  normal¬ 
ized  form, 

f  P(z,6.<t>)d(l  -  1  (25) 

4ir 

The  computational  scheme  offers  several  options  for 
specification  of  the  average  single  scattering  phase  func¬ 
tion  P(zfi),  for  each  designated  atmospheric  layer. 
Detailed  estimates  of  PM(fi)  available  from  prescribed 
aerosol  models  such  as  those  associated  with  LOWTRAN 
(Shettle  and  Fenn,  1979),  may  be  entered  directly  and 
used  in  a  table  look-up  format.  A  second  option  is  the 
representation  of  the  aerosol  phase  functions  by  two-term 
Henyey-Greenstein  functions  (Irvine,  1968)  as  follows, 

P» 03,y,,^2,c)  -  c/>//g(0,£|)  +  (1  -  c)PMc(j 3,g:)  ,  (26) 

where 

PHG(0.g)  -  (1  -  g2)/  1 4ir (1  -  2#cos')  +  s2)VJl  .  (27) 
and  the  asymmetry  factor,  g ,  is  given  by 


g  -  4  f  f><0)cos/3  sin))  d  0  (28) 

2  o 

To  the  extent  that  the  Henyey-Grdenstein  asymmetry  fac¬ 
tors,  *|  and  g 2,  and  the  partitioning  factor,  c.  can  be 
estimated  from  reference  atmosphere  calculations  they  can 
be  entered  for  the  specification  of  the  required  scattering 
phase  functions.  In  the  absence  of  information  that  may 
serve  to  identify  the  directional  scattering  properties  of 
individual  aerosol  layers,  yet  another  option  may  be  used 
as  described  in  Hering  (1981).  Empirical  functions  were 
developed  which  prescribe  the  Henyey-Greenstein  func¬ 
tion  parameters,  g)t  g2,  and  c,  as  a  function  of  the  scatter¬ 
ing  ratio  Q(z).  The  derived  expressions  are  based  to  a 
large  extent  on  the  average  phase  functions  for  selected 
ranges  of  scattering  coefficient  as  measured  by  Barteneva 
(1960) 

The  other  direct-entry  variable  for  each  designated 
atmospheric  layer  is  the  single  scattering  albedo,  o>(z),  It 
is  defined  by  the  expression 


4.0  APPLICATION  TO  OVERCAST  SKY  CONDITION 

The  clear-sky  modeling  concepts  and  computational 
scheme  can  be  applied  in  the  same  form  to  calculate  the 
visible  spectrum  radiance  fields  and  contrast  transmittance 
for  continuous  cloud  cover  conditions  Given  reliable 
empirical  estimates  of  cloud  optical  properties,  the  delta- 
Eddington  method  yields  good  estimates  of  the  irradiance 
profiles  in  overcast  atmospheres  In  turn,  the  directional 
path  radiance  fields  and  the  contrast  transmittance  for 
selected  slant  paths  may  be  calculated  from  the  expres¬ 
sions  given  in  Section  2  In  practice,  one  must  rely  on 
cloud  model  estimates  of  the  optical  depth,  the  single 
scattering  phase  function,  and  the  single  scattering  albedo 
as  a  function  of  standard  cloud  type. 


4.1  Estimates  of  Optical  Depth  from 
Conventional  Cloud  Observations 

Rough  empirical  approximations  of  the  average  opti¬ 
cal  depth  as  a  function  of  cloud  type  were  derived  through 
reference  to  a  climatological  summary  of  the  SOLMET 
(1977)  data  base  as  included  in  a  report  by  Shapiro 
(1982).  The  SOLMET  data  base  was  compiled  by  the 
National  Climatic  Center  under  the  sponsorship  of  the 
Department  of  Energy.  It  consists  of  a  long  record  of 
edited  and  adjusted  solar  radiation  data  (hourly  averages) 
for  26  National  Weather  Service  stations  spread 
throughout  the  United  States  coupled  with  the  conven¬ 
tional  hourly  surface  observations  for  these  stations.  The 
climatological  tabulations  of  the  broad  band  solar  irradi¬ 
ance  were  available  in  the  form  of  averages  for  all  stations 
and  all  seasons  for  standard  cloud  types  (10/10  cloud 
cover)  for  selected  solar  zenith  angles  These  data  were 
combined  with  the  average  irradiance  data  for  cloudless 
skies  to  calculate  the  ratio  of  the  average  broad  band  irra¬ 
diance  for  the  sky  covered  by  clouds  of  specific  types  to 
the  average  broad-band  irradiance  for  cloudless  atmo¬ 
spheres. 

It  is  important  to  note  that  the  calculated 
overcast/clear  irradiance  ratios  for  different  cloud  types 
and  solar  angles  derived  from  the  SOLMET  data  base 
agree  very  closely  with  the  average  overcast/clear  irradi¬ 
ance  ratios  calculated  by  (Haurwitz  (1948))  for  compar¬ 
able  cloud  classifications  using  8  years  of  observations  at 
the  Blue  Hill  Observatory  in  Massachusetts. 


w(z )  “  s  (z)/ a  (z  )  . 


(29) 


To  the  extent  the  fractional  transmittance  for  the 
different  cloud  types  as  determined  from  the  broed-band 


(. 295-2. IVtm)  pyronometer  measurements  are  representa¬ 
tive  in  general  of  the  visible  spectrum,  one  can  use  the 
delta- Eddington  method  to  estimate  the  effective  average 
optical  depth  associated  with  the  individual  cloud  types. 
The  results  of  trial  calculations  made  to  obtain  rough  esti¬ 
mates  of  the  optical  depth  of  standard  cloud  types  are 
shown  in  Fig.  6.  Through  successive  approximation, 
determination  was  made  of  the  best  least-squares  fit 
between  the  average  overcast/clear  irradiance  ratios  as  a 
function  of  solar  zenith  angle  derived  from  the  SOLMET 
data  base  and  the  corresponding  values  calculated  with  the 
delta-Eddington  model.  A  central  wavelength  of  550  nm 
and  an  average  surface  reflectance  of  .07  were  used  for  the 
model  calculations  and  absorption  within  the  cloud  layer 
was  assumed  negligible.  The  altitude  limits  of  the  cloud 
layers  were  chosen  so  as  to  conform  in  general  with  typical 
values  for  each  cloud  type. 

Except  as  noted  below,  the  asymmetry  factor,  g, 
was  assumed  equal  to  .85  in  calculations  for  all  cloud 
types.  It  should  be  emphasized  that  the  model  calcula¬ 
tions  are  sensitive  to  variations  in  g.  However,  in  the 
case  of  water  clouds,  the  variation  in  g  is  small  (Twomey, 
1980)  and  the  resultant  uncertainties  are  negligible  in 
comparison  with  the  natural  variations  in  optical  depth  for 
each  cloud  type,  fn  the  case  of  ice  clouds,  very  significant 
fluctuations  in  g  are  associated  with  variations  in  ice  cry¬ 
stal  size,  shape,  and  orientation.  For  this  reason,  as 
emphasized  by  Welch,  el  al.  1980,  attempts  to  estimate 
the  optical  depth  of  clouds  composed  of  non-spherical  par¬ 
ticles  from  observed  cloud  reflectance  (or  transmittance) 
are  subject  to  considerable  error.  For  comparison  pur¬ 
poses,  estimates  were  derived  for  both  a  nominal  value  of 
.85  for  the  asymmetry  factor  in  cirrus  clouds  and  a  marked 
reduction  to  .75.  In  agreement  with  the  results  of  Welch, 
et  al.  1980,  the  calculated  optical  depths  corresponding  to 
the  same  average  cirrus  cloud  transmittance  differ  by  more 


ns-  f.  Avenge  overcan/deer  irradiance  ntios  from  SOLMET 
data  twee  (aotid  lines)  end  estimated  with  delta-Eddington  approx¬ 
imation  (dashed  lines).  Listed  are  the  values  of  optical  depth,  r, 
living  the  best  correspondence. 


Table  3. 

The  approximate  average  and  ±  standard  deviation 
values  of  optical  thickness  for  various  cloud  forms  as 
calculated  using  the  delta- Eddington  approximation  and 
the  SOLMET  data  base  The  assumed  value  of  asym¬ 
metry  lactor  is  given  in  column  g  For  reference  pur¬ 
poses  a  value  of  scattering  ratio,  Q(550nm),  is  shown 
for  each  value  of  dr  as  determined  for  a  cloud  of 
geometric  cloud  thickness  of  1  km  and  a  base  altitude  as 
indicated. 


Cloud  Type 

- r 

Optical  I  QI550  nm) 
Thickness  j  -  1  km 

Base  | 

All  (km)  |  g 

1 

Sl/Sc 

22  0 

2200 

1  km 

13  5 

1350 

85 

410 

4100 

Cu/Cb 

150 

1500 

1  km 

100 

1000 

85 

26  5 

2650 

As/  Ac 

IS 

1950 

3  km 

98 

1200 

85 

25  3 

3100 

Ns  (precip.) 

43  0 

4300 

1  km 

1 

207 

2700 

85 

126 

12600 

Ci/Cs  (thick) 

44 

750 

6  km 

i 

32 

550 

.85 

62 

1050 

j 

Ci/Cs  Uhin) 

15 

250 

6  km 

07 

125 

85  1 

24 

400 

• 

Ci/Cs  (thick) 

26 

440 

6  km 

15 

250 

.75 

43 

725 

Ci/Cs  (thin) 

O' 

120 

6  km 

02 

30 

75 

1  5 

250 

than  a  factor  of  2  for  the  two  assumed  values  of  g  as 
shown  in  Table  3. 

As  indicated  by  the  examples  in  Fig.  6,  one  can 
obtain  a  rather  close  fit  of  the  average  irradiance  ratios  vs 
solar  zenith  angle  for  the  standard  cloud  types  using  the 
delta-Eddington  model.  The  optical  depths  corresponding 
to  the  closest  approximations  of  the  SOLMET  data  aver¬ 
ages  for  each  cloud  type  are  listed  in  Table  2.  Rough  esti¬ 
mates  of  the  variability  in  the  optical  depth  were  obtained 
in  a  similar  manner  from  the  SOLMET  data  The  tabu¬ 
lated  values  of  optical  depth  labeled  "±  rms  variation’ 
were  calculated  from  SOLMET  irradiance  ratios  equal  to 
the  average  ratio  for  a  given  cloud  type  and  solar  zenith 
angle  plus  and  minus  the  root-mean-square  departure  of 
the  individual  SOLMET  values  of  irradiance  ratio  from  the 
averages  for  that  cloud  type. 
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Thus,  for  overcast  sky  conditions,  the  radiance  field 
at  any  altitude  may  be  calculated  using  the  optical  proper¬ 
ties  given  in  Table  2  corresponding  to  the  observed  cloud 
type  and  altitude  and  the  estimated  input  variables  for  the 
additional  atmospheric  layers  above  and  below  the  cloud 
layer.  The  optical  scattering  ratio  for  the  cloud  layer  is 
given  by  the  cloud  optical  depth  divided  by  the  Rayleigh 
optical  depth  for  the  corresponding  altitude  interval  of  the 
cloud.  In  view  of  the  uncertainty  in  the  asymmetry  factor 
for  cirrus  cloud  particle  distributions,  estimates  for  g  equal 
to  75  are  shown  as  well  as  for  the  more  representative 
value  of  .85.  In  the  absence  of  other  information,  opera¬ 
tional  use  of  the  results  for  g«.85  are  suggested 

Several  comparisons  of  the  calculated  spectral  radi¬ 
ance  fields  with  the  high  resolution  aircraft  measurements 
for  both  clear  and  cloudy  atmospheres  are  given  in  the  fol¬ 
lowing  section. 

5.0  COMPARISONS  OF  MODEL  RESULTS 
WITH  MEASURED  DATA 

From  the  fundamental  Eq.  4,  we  note  that  the  relia¬ 
bility  of  model  calculations  of  slant-path  contrast  transmit¬ 
tance  are  dependent  directly  upon  the  accuracy  with  which 
the  radiance  fields  L,  (z,0,<£)  are  determined.  For  this 
reason,  continuing  emphasis  has  been  given  to  direct  com¬ 
parisons  of  the  radiance  fields  as  calculated  by  the  model¬ 
ing  procedures  and  the  corresponding  radiance  fields 
observed  as  part  of  the  extensive  airborne  measurement 
program.  The  upper  and  lower  scanning  radiometers 
measured  the  apparent  radiance  fields  surrounding  the  air¬ 
craft  as  it  probed  the  lower  troposphere  up  to  an  altitude 
of  6  km.  The  processed  radiometer  measurements  pro¬ 
vided  data  with  5-degree  resolution  in  the  4  specific 
wavebands  commensurate  with  the  integrating  nephelome- 
ter  measurements  of  total  volume  scattering  coefficient 
(peak  wavelengths  of  475,  550,  670  and  750  nm). 

Several  of  the  comparisons  between  model  calcula¬ 
tions  and  measured  radiance  fields  are  shown  in  Figs.  8  to 
12.  The  plotted  data  are  presented  for  either  one  of  two 
vertical  planes  (0-180  deg  or  90-270  deg)  with  respect  to 
the  solar  azimuth  and  depict  data  sweeping  from  the  nadir 
(9  —  180  deg)  up  through  the  horizon  and  zenith  and 
then  downward  through  the  opposite  horizon  to  nadir,  as 
shown  in  Fig.  7.  A  series  of  examples  were  selected  in 
order  to  illustrate  model  performance  over  a  wide  range  of 
haze-layer  extinction,  solar  zenith  angle,  and  surface 
reflectance  (land  and  sea)  for  both  clear  and  overcast  sky 
conditions. 

The  specific  input  data  for  each  of  the  illustrated 
examples  are  listed  in  Appendix  B.  Except  for  the  scatter¬ 
ing  mixing  ratio  profile  in  the  lower  troposphere,  the  other 
input  variables  were  approximated  by  their  average  or 
climatological  values.  The  scattering  mixing  ratio  for  the 
boundary  layer  was  entered  as  determined  directly  from 
the  nephelometer  measurements  for  that  particular  flight 
profile.  The  scattering  ratio  for  the  upper  troposphere  and 
stratosphere  during  this  period  of  minimal  volcanic  activity 
was  held  constant  for  all  cases  and  set  equal  to  a  charac¬ 
teristic  value  of  1.2  at  550  nm.  The  mixing  ratio  for  other 


Fl*.  7.  Coordinate  system  for  the  calculation  of  sky  and  terrain 
radiance  and  slant-path  contrast  transmittance. 


wavelengths  at  upper  levels  was  determined  assuming  the 
wavelength  dependence  of  the  aerosol  scattering 
coefficient  equal  to  (\/X550)  1  The  aerosol  single  scatter¬ 
ing  phase  function  for  all  haze  layers  and  all  cases  was 
specified  by  a  two-term  Henyey-Greenstein  function  with 
parameters;  £,  -  714,  g }  -  -.613,  and  r  -  .963  The 
corresponding  values  for  the  cloud  layers  are  .890,  -.660 
and  .980. 

Optical  depths  for  the  overcast  sky  layers  were 
specified  equal  to  the  average  value  for  the  observed  cloud 
type  as  developed  in  Section  4  (see  Table  2).  The  single 
scattering  albedo  and  surface  reflectance  were  assumed 
equal  to  average  values  for  the  prevailing  conditions  For 
flights  conducted  over  open  water,  the  average  Fresnel 
reflectance  distribution  for  wind-ruffled  sea  (see  Sec¬ 
tion  3.1)  was  used. 


5.1  Comparative  Results  for 
Individual  Flights 

Right  C-379,  Fig  8,  was  a  midday  flight  ofT  the 
south  coast  of  Lolland  Island,  Denmark.  The  example 
depicts  the  characteristic  changes  with  altitude  of  the 
upper  and  lower  hemisphere  radiance  fields  under  cloud¬ 
less  skies  and  light  haze  conuitions.  Note  the  characteris¬ 
tic  uniformity  and  symmetry  of  the  sea  surface  radiance 
fields  in  the  cross-sun  plane. 

Right  C-466,  Fig.  9,  was  conducted  over  open  farm¬ 
land  and  wooded  areas  in  northwestern  Germany.  Scat¬ 
tered  cirrus  clouds  and  isolated  patches  of  altocumulus 
were  observed.  A  layer  of  moderate  haze  extended  to  an 
altitude  of  1350m.  Notice  the  marked  fluctuations  in  thr 
observed  apparent  terrain  radiance  associated  with  vari 
tions  in  ground  cover  (green  and  brown  fields  and  dark 
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Flf.  9.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  1981d)  (solid  line)  and  model  cal 

culations  (dashed  line). 


woods),  as  opposed  to  the  smooth  terrain  radiance  distri¬ 
bution  calculated  with  the  assumption  of  uniform  Lamber¬ 
tian  surface  reflectance. 

Right  C-469,  Figs.  10  and  11,  was  an  afternoon 
flight  over  open  farmland  in  central  Netherlands.  The 
examples  illustrate  the  variations  in  radiance  fields  with 
altitude  in  dense  haze  conditions.  Depicted  also  are  the 
characteristic  differences  in  the  radiance  distributions  as 
observed  in  the  narrow  blue  passband  (475  nm)  as  shown 
in  Fig.  10  and  in  the  near  infrared  passband  (750  nm) 
shown  in  Fig.  1 1 .  Note  in  particular  the  reduction  in  over¬ 
head  sky  radiance  at  750  nm  associated  with  diminished 
molecular  and  aerosol  particle  scattering  and  the  enhanced 
apparent  terrain  radiance  due  to  the  greater  surface 
reflectance  in  the  near  infrared  portion  of  the  spectrum. 

Right  C-401,  Fig.  12,  was  a  winter  day  flight  with 
low  sun  elevation  over  a  rural  area  in  northwestern 
France.  Scattered  cumulus  and  scattered  thin  cirrus  clouds 
were  observed  with  a  layer  of  light  haze  extending  to  an 
altitude  of  about  1  km.  The  relatively  bright  horizon  in 
contrast  with  low  overhead  sky  and  terrain  radiance  is  a 
prominent  feature  of  low-sun  radiance  fields.  As  indicated 
by  this  example,  the  good  agreement  between  model  cal¬ 
culations  and  measured  radiance  fields  characteristically 
found  in  the  comparisons  extends  to  the  cases  with  low 
solar  elevation  angle. 

Flight  C-422,  Fig.  13,  is  also  an  over  water  flight  off 
the  coast  of  Denmark,  but  on  this  day  the  sky  was  over¬ 
cast  with  stratocumulus  clouds  having  a  base  altitude  of 
1800m.  Heavy  haze  was  present  at  all  levels  below  the 
cloud  base.  The  example  illustrates  the  highly  symmetric 
and  uniform  distribution  of  cloud  and  sea  radiance  typical 
of  measurements  made  with  the  flight  path  over  open 
water  and  under  an  overcast  cloud  layer. 

Flight  C-465,  Fig.  14,  was  an  afternoon  flight  in  mid 
summer  over  northwestern  Germany.  The  sky  was  over¬ 
cast  with  stratocumulus  clouds  having  a  base  altitude  of 
1600m.  A  layer  of  light  haze  was  observed  below  the 
cloud  base.  In  contrast  to  the  previous  case,  the  radiance 
field  was  more  irregular  in  structure  due  to  the  variations 
in  cloud  optical  thickness  and  to  the  variations  in  ground 
cover  reflectivity. 

Right  C-435,  Fig.  15,  was  a  winter  flight  over  a 
snow  covered  plateau  in  southwestern  Germany.  The  sky 
was  overcast  with  altostratus  clouds  with  a  base  altitude  of 
3  km.  The  apparent  radiance  of  the  snow  surface  as 
measured  at  an  altitude  of  6  km  is  roughly  equal  to  the 
apparent  radiance  of  the  overcast  sky.  The  minimum 
values  of  lower  hemisphere  radiance  are  associated  with 
underlying  areas  of  dark  woods  having  an  effective  aver¬ 
age  reflectance  of  about  .04  in  this  case.  Notice  that  the 
effective  reflectance  of  the  snow  covered  areas  is  some¬ 
what  higher  than  the  average  value  of  0.75  that  was 
assumed  for  the  model  calculations. 

The  correspondence  between  the  field  measure¬ 
ments  and  the  model  calculations  illustrated  above  are 
representative  in  general  of  the  many  such  comparisons 
made  during  the  course  of  model  development  and  valida¬ 


tion.  It  should  be  emphasized  that  the  accuracy  of  the 
illustrated  model  calculations  is  enhanced  considerably  by 
the  availability  of  the  detailed  vertical  profile  of  scattering 
coefficient  in  the  boundary  layer,  which  is  the  msgor  deter¬ 
minant  of  the  fluctuations  in  path  radiance  in  the  tropo¬ 
sphere.  On  the  other  hand,  the  correspondence  between 
the  actual  and  calculated  radiance  fields  could  be  improved 
through  relaxation  of  the  constraint  that  other  input  vari¬ 
ables  be  held  equal  to  their  typical  or  average  values.  For 
example,  identification  of  an  aerosol  model  appropriate  for 
the  boundary  layer  in  a  given  synoptic  situation  would  pro¬ 
vide  estimates  of  the  single  scattering  phase  function  and 
single  scattering  albedo  that  are  more  realistic  than  the 
prescribed  overall  average  values  used  in  illustrated  com¬ 
parisons.  A  brief  assessment  of  the  sensitivity  of  esti¬ 
mates  of  the  maximum  detection  range  of  distant  objects 
to  uncertainties  in  both  atmospheric  and  target  factors  is 
given  in  the  following  Section. 


6.0  SENSITIVITY  OF  TARGET  DETECTION  RANGE 
TO  ENVIRONMENTAL  AND  TARGET  FACTORS 

Diagnostic  calculations  have  been  made  routinely 
during  the  course  of  the  investigations  to  determine  the 
sensitivity  of  the  model  calculations  to  variations  in  the 
input  data.  In  prior  investigations  (Hering,  1981),  atten¬ 
tion  was  directed  primarily  toward  analysis  of  the  sensi¬ 
tivity  of  contrast  transmittance  to  changes  or  uncertainties 
in  the  relevant  optical  properties  of  the  atmosphere.  The 
interpretations  of  these  results  are  straightforward  in  that, 
except  for  scintillation  effects,  Eq.  4  describes  completely 
the  atmospheric  effects  upon  contrast  transmittance.  On 
the  other  hand,  contrast  transmittance  must  be  coupled 
with  the  effects  of  specific  target  factors  to  estimate  the 
maximum  distance  an  object  can  be  detected  and 
identified.  A  given  change  in  the  slant-path,  contrast 
transmittance  may  or  may  not  be  critical  for  target  detec¬ 
tion  depending  the  minimum  contrast  required  for  detec¬ 
tion.  Thus,  the  assessment  of  the  combined  influences  of 
atmospheric  and  target  factors  is  important. 

For  this  purpose,  analytic  techniques  were  intro¬ 
duced  into  the  computer  program  to  calculate  visual  detec¬ 
tion  range  as  an  additional  output  variable.  The  represen¬ 
tation  of  target  information  and  visual  search  factors  in 
this  preliminary  analysis  are  greatly  simplified.  Neverthe¬ 
less,  the  results  provide  important  insight  into  the  sensi¬ 
tivity  of  visual  detection  distance  to  changes  in  specific 
atmospheric  variables  for  selected  combinations  of  target 
characteristics. 

6.1  Determination  of  Visual 
Detection  Range 

The  analytic  representations  of  the  contrast  thres¬ 
holds  for  visual  detection  were  derived  from  basic  visual 
data  from  the  Tiffany  experiments  (Blackwell,  1946)  and 
the  Visibility  Laboratory  experiments  (Taylor,  1 964)  which 
were  adapted  and  summarized  by  Gordon  (1979).  The 
graphical  representations  of  threshold  contrast  as  a  func¬ 
tion  of  angular  target  size  shown  in  Fig.  16  are  taken 
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Fig.  10.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  I981d)  (solid  line)  and  model  cal' 

culations  (dashed  line) 
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rig.  11.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  1981d)  (solid  line)  and  model  cal 

culations  (dashed  line). 
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Fig.  12.  C  omparison  of  sky  and  lerrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  1 98 Id)  (solid  line!  and  model  cal¬ 
culations  (dashed  line) 


Fig.  13.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  198  Id)  (solid  line)  and  model  cal¬ 
culations  (dashed  line). 
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Fig.  14.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  W81dl  t solid  line •  and  mode!  cal¬ 
culations  (dashed  line ) 


C-435  BIRKHOF.  W  GERMANY 
ALTOSTRATUS  OVERCAST 
SNOW  COVERED  PLATEAU 
550nm 
-  65  deg 


Fig.  15.  Comparison  of  sky  and  terrain  radiance  as  measured  by  airborne  scanning  radiometer  (Johnson  1 98 1  d )  (solid  line)  and  model  cal¬ 
culations  (dashed  line) 
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Fig.  16.  Angular  subtense  of  larger  as  a  function  of  threshold 
contrast  for  a  99  percent  probability  (confidence)  of  detection  and 
a  lack  of  knowledge  of  target  position  of  ±  4  deg  or  more  (Gor¬ 
don,  1979)  bB,  is  the  adaptation  brightness  and  Si  is  the 
assumed  visual  fixation  (dwell)  time 


directly  from  Fig  1  of  Gordon  (1979).  Two  search  modes 
are  represented  The  first  corresponds  to  long  duration 
viewing  0  24  seconds.  Tiffany  data).  The  second 
corresponds  to  a  viewing  time  of  one-third  of  a  second 
(Visibility  Laboratory  data),  which  is  the  approximate 
fixation  time  associated  with  search  performed  by  the 
human  eye  (White,  1964). 

From  Fig  16  we  see  that,  for  the  normal  range  of 
daylight  conditions,  the  visual  contrast  threshold  does  not 
vary  significantly  with  background  luminance  However, 
the  contrast  threshold  does  vary  markedly  with  the  angular 
subtense  of  the  target.  The  relationships  shown  in 
Fig.  16,  represent  the  contrast  threshold  associated  with 
confident  detection  (99  percent  probability)  under  condi¬ 
tions  where  the  lack  of  knowledge  of  target  position  is  ±  4 
degrees  or  more. 

For  application  to  the  problem  at  hand,  the 
prescribed  relationships  given  in  Fig.  16  were  approxi¬ 
mated  analytically  by  the  expression. 


lr  y  -  A  +  B  (-In  Cr)  +  C(-ln  CT)2  +  C(-ln  Cr>5  (30) 


of  the  model  calculations  of  slant  path  contrast  transmit¬ 
tance.  The  computational  scheme  provides  discrete  values 
of  contrast  transmittance  for  designated  paths  of  sight.  As 
an  additional  option,  the  contrast  transmittance  of  air  to 
ground  observation  paths  for  fixed  values  of  azimuthal 
viewing  angle  and  sensor  altitude  is  represented  through 
the  method  of  least  squares  as  a  continuous  function  of 
path  length  (range)  by  the  expression, 

C  C 

r  -  a  +  *|ln  -pr-  (r,fl)l  +  c (In  (z,0)P  (31) 


Since  the  angular  subtense  of  the  target  is  less  than  a 
degree  or  so.  the  slant  range,  r,  can  be  approximated  by 
r  =  d/y  where  d  is  the  effective  diameter  of  the  target 
Given  the  inherent  target  contrast,  C„,  the  range  associ¬ 
ated  with  a  selected  value  of  apparent  contrast  O.  is 
determined  by  Eq  31.  The  maximum  detection  range  for 
specified  sensor  altitude  and  azimuthal  viewing  angle 
corresponds  to  the  range  where  C,  =  C, .  The  method  of 
successive  approximation  yields  rapidly  converging  solu¬ 
tions  for  maximum  detection  range  from  trial  values  of  C, 
and  C,  using  Eqs  30  and  31 

6.2  Diagnostic  Calculations  of 
Visual  Detection  Distance 

A  series  of  additional  model  calculations  were  car¬ 
ried  out  to  illustrate  the  relative  sensitivity  of  target  detec¬ 
tion  range  to  changes  or  uncertainties  in  some  of  the  key 
variables.  For  this  purpose,  reference  conditions  were 
assumed  which  reflect  a  typical  visual  detection  scenario 
for  hazy  atmospheres  in  the  daytime.  In  an  attempt  to 
isolate  the  effects  of  fluctuations  in  individual  factors,  the 
procedure  was  to  vary  a  given  parameter  over  its  normal 
range  while  holding  all  other  reference  conditions  con¬ 
stant.  The  reference  atmospheres  used  in  this  analysis  are 
summarized  in  Appendix  B  (Ref.  A  for  clear  sky  and 
Ref.  B  for  overcast  sky).  The  reference  atmospheres  have 
a  surface  equivalent  visibility  of  20  km  and  a  boundary- 
layer  hazy  depth  of  1  km.  Except  as  noted,  an  ideal  circu¬ 
lar  target  is  assumed  having  a  diameter  of  30m  and  an 
inherent  contrast  of  one.  The  functional  representation  of 
psychophysical  data  of  Taylor  (1964)  corresponding  to  the 
0.33  sec  fixation  time  (see  Fig.  16)  was  used  for  the  trial 
calculations  except  for  the  one  comparative  example  (para¬ 
graph  6.2.5)  based  upon  the  Tiffany  visual  threshold  data 
for  a  fixation  time  greater  than  24  seconds. 

Visual  Horizontal  Detection  Distance 
versus  Observer  Altitude  and  Azimuthal 
Viewing  Angle  (Fig.  17) 


where  y  is  the  angular  size  of  the  target.  C>  is  the  con¬ 
trast  detection  threshold  and  A,B,C  and  D  are  the  parame¬ 
terization  coefficients  as  determined  for  the  two  search 
modes 

Given  the  above  method  for  the  representation  of 
target  factors,  let  us  now  consider  analytic  representation 


For  intermediate  size  targets  of  good  inherent  con¬ 
trast,  the  maximum  horizontal  detection  distance  (ground 
range)  increases  significantly  with  increasing  observer  alti¬ 
tude  in  the  lower  troposphere  above  the  toundary  layer. 
The  variations  in  ground  range  with  azimuthal  viewing 
angle  are  determined  largely  by  the  directional  dependence 
of  the  path  radiance  as  defined  by  the  combined  phase 
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«*  17.  Visual  horizontal  detection  distance  as  a  function  of  Fig.  18.  Visual  horizontal  detection  distance  as  a  function  of 
observer  altitude  and  azimuthal  observation  angle  inherent  contrast  for  30-m  target 


function  for  molecular  and  aerosol  scattering.  For  the  typ¬ 
ical  condition  in  the  free  atmosphere  above  the  haze  layer 
where  the  aerosol  particles  have  low  concentration,  the 
difference  in  detection  distance  in  the  upsun  (<t>  —  0  deg) 
and  downsun  (d>  —  180  deg)  viewing  directions  remains 
rather  constant  with  increasing  observer  altitude  above  the 
boundary  layer  and  increases  with  increasing  solar  zenith 
angle.  Of  special  interest  in  Fig.  17  is  the  relatively  small 
variation  in  detection  range  with  azimuthal  angle  in  the 
antisolar  hemisphere  as  indicated  by  the  rather  close 
correspondence  in  the  180  degree  (downsun)  and  90 
degree  (cross  sun)  profiles.  It  should  be  noted  that  the 
trial  calculations  which  follow  in  subsequent  examples 
assume  a  90-degree  azimuthal  viewing  direction. 

Visual  Detection  Distance  versus 
Inherent  Contrast  (Fig.  18). 

The  rate  of  increase  in  visual  detection  range  with 
increasing  observer  altitude  depends  significantly  upon  the 
initial  target  contrast.  Notice  in  Fig.  18  the  dramatic 
increase  in  detection  range  with  increasing  viewing  altitude 
for  bright  targets  viewed  against  a  dark  background 
(C„  -  5)  as  contrasted  with  the  restricted  increase  in 
detection  distance  for  low  contrast  targets  (C„  -  0.2)  As 
a  point  of  reference  information,  the  detection  range  is 
independent  of  the  sign  of  the  target  contrast.  For  targets 
darker  than  the  background,  the  lower  contrast  limit  is 
C»  —  - 1  (black  target) .  The  inherent  contrast  of  targets 
brighter  than  background  is  unbounded.  The  assumed 
reference  value  of  target  contrast  of  C0  -  1  (dashed 
curve)  corresponds  to  the  case  where,  under  the  same 
irradiance,  the  reflectivity  of  the  target  is  twice  the  back¬ 
ground  reflectivity. 

Visual  Detection  Distance  versus 
Inherent  Contrast  for  Reduced 
Target  Size  (Fta.  19). 

Diagnostic  model  calculations  for  a  smaller  5-m  tar¬ 
get  diameter  (Fig.  19)  show  the  corresponding  reduction 


Fig.  19.  Visual  horizontal  detection  distance  as  a  function  of 
inherent  contrast  for  5-m  target 

in  visual  detection  distance  as  compared  with  the  results 
for  the  30-m  target  (Fig.  18).  Notice,  for  example,  in  the 
case  of  a  target  of  0.2  inherent  contrast  under  the  assumed 
conditions,  there  is  little  advantage  in  increased  detection 
range  to  be  gained  by  climbing  to  an  observation  altitude 
above  the  haze  layer  top  and  that  the  reference  target  hav¬ 
ing  0,2  contrast  cannot  be  detected  with  99  percent  proba¬ 
bility  above  a  viewing  altitude  of  6  km.  The  calculated 
visual  detection  range  for  5-m  targets  for  this  example  is 
everywhere  less  than  12  km  for  targets  having  an  inherent 
contrast  less  than  one. 

Visual  Detection  Distance  as 
a  Function  of  Environmental 
Factors  (Fig.  20). 

Sensitivity  calculations  of  the  response  of  ground 
detection  range  to  selected  changes  in  individual  factors 
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Fig.  20.  Sensitivity  of  horizontal  detection  distance  to  changes  in  relevant  factors.  Solid  bar  in  each  case  (22  km)  corresponds  to  the  calcu¬ 
lation  for  reference  target  (30-m  diameter)  and  reference  atmosphere  parameters.  Surface  target  with  a  contrast  of  one  is  viewed  against  at 
a  uniform  background  from  an  altitude  of  3  km  For  each  case,  all  variables  are  fixed  equal  to  reference  values  except  for  the  variable 
under  examination 


are  shown  in  Fig.  20  for  a  ground  target  viewed  from  an 
altitude  of  3  km.  All  variables  except  the  factor  under 
examination  are  held  constant  and  equal  to  the  reference 
values  listed  in  paragraph  6.1  and  Appendix  B  (Ref.  A). 
The  scattering  and  absorption  coefficients  are  allowed  to 
change  only  within  the  primary  haze  layer.  The  calcula¬ 
tions  assume  a  A“'  wavelength  dependence  of  the  aerosol 
scattering  coefficient  and  no  changes  in  surface  reflectance 
with  wavelength.  The  illustrated  changes  in  detection  dis¬ 
tance  relative  to  changes  in  single  scattering  albedo  were 
calculated  assuming  that  the  optical  thickness  of  the  haze 
layer  remains  constant,  i.e.  the  scattering  ratio  was 
changed  by  an  amount  sufficient  to  compensate  for  the 
postulated  change  in  absorption  coefficient 

The  solid  bar  in  each  of  the  graphs  shown  in  Fig.  20 
depicts  the  ground  range  (22  km)  associated  with  the 
reference  values  for  the  individual  variables.  Notice  for 
the  prescribed  conditions  that  a  decrease  in  target  distance 
of  about  25  percent  is  associated  with:  (Da  reduction  in 
target  size  by  50%,  (2)  a  reduction  in  surface  equivalent 
visibility  by  35%;  (3)  an  increase  in  solar  zenith  angle 
from  60  to  78  degrees.  (4)  a  decrease  in  surface 
reflectance  by  50  percent,  and  (5)  a  change  in  wavelength 
from  550  to  450  nm.  The  detection  distance  is  less  sensi¬ 
tive  to  changes  in  the  absorption  contribution  to  a  fixed 
value  of  total  volume  extinction  coefficient  in  the  boun¬ 
dary  layer  It  should  be  emphasized  that  the  results 
shown  in  Fig.  20  are  specific  only  for  the  conditions 
assumed  for  the  diagnostic  calculations.  Simultaneous 
changes  in  two  or  more  variables  could  of  course  result  in 


a  much  broader  range  in  the  variation  of  target  detection 
distance. 

Visual  Detection  Distance  versus 
Visual  Fixation  Time  (Fig.  21). 

The  variations  in  target  detection  distance  that  are 
associated  with  changes  in  stimulus  duration  are  shown  in 
Fig.  21.  The  two  functional  relationships  for  contrast 
threshold  data  used  for  the  calculations  are  the  same  as 
those  depicted  in  Fig.  16.  As  discussed  above,  the  first 
applies  to  unlimited  viewing  time  (>  24  sec)  and  the 
second  applies  to  the  average  fixation  time  (1/3  sec)  for 
the  performance  of  visual  search  by  the  human  eye. 
Comparative  calculations  for  the  two  viewing  conditions 
are  shown  for  two  atmospheres.  Both  atmospheres  and 
target  factors  are  the  same  as  reference  conditions  cited 
above  except  that  the  haze  layer  of  depth  1.3  km  is 
represented  by  the  LOWTRAN  (Shettle  and  Fenn,  1979) 
rural  model  atmosphere  for  relative  humidity  equal  to  70 
percent  for  the  one  case  and  for  relative  humidity  equal  to 
99%  in  the  other  case.  Notice  that  the  relative  difference 
in  detection  distance  versus  glimpse  time  for  the  example 
with  high  humidity  and  high  extinction  coefficient  is 
roughly  the  same  as  that  for  the  lower  humidity  and  lower 
extinction  coefficient  example. 

Visual  Detection  Distance 
for  Overcast  Sky  (Fig.  22). 

A  cloudless  atmosphere  was  assumed  for  the  sensi¬ 
tivity  analyses  described  above.  Let  us  now  consider 
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n«.  21.  Visual  horizontal  detection  distance  versus  altitude  for 
fixation  times  (At)  of  0.33  sec  and  infinity.  Calculations  assume 
LOWTRAN  rural  model  for  boundary  layer  with  70  percent  humi¬ 
dity  (equivalent  surface  visibility  of  about  20  km)  and  with  99  per¬ 
cent  humidity  (equivalent  surface  visibility  of  about  5  km) 


Fig.  22.  Comparison  of  calculated  target  detection  distance  for  a 
reference  atmosphere  with  and  without  an  overcast  cloud  layer 
Calculations  for  clear  sky  case  (solid  lines)  are  shown  for  solar 
zenith  angles  of  25,  65,  and  85  deg.  Calculations  for  case  with 
stratocumulus  overcast  (dashed  line)  are  for  solar  zenith  angle  of 
55  deg. 


briefly  the  behavior  of  contrast  transmittance  and  the  vari¬ 
ations  in  target  detection  distance  that  are  associated  with 
overcast  sky  conditions.  In  the  region  below  the  cloud 
layer  in  the  daytime,  both  the  inherent  background  radi¬ 
ance  „/*(/, ,0,<6)  and  the  apparent  radiance  bL,U  ,9,  <t>) 
are  essentially  proportional  to  the  downwelling  radiance  at 
the  cloud  base  when  the  direct  sun  component  of  the  irra- 
diance  is  small  in  comparison  with  the  diffuse  component. 
Thus,  except  for  thin  clouds  with  significant  direct  sun 
penetration,  the  contrast  transmittance  of  slant  paths 
below  the  cloud  deck  is  to  a  good  approximation  indepen¬ 
dent  of  the  downwelling  irradiance  from  the  cloud  layer 
(see  Eq.  4).  However,  the  path  radiance  and  radiance 
transmittance  (hence  the  contrast  transmittance)  are  of 
course  strongly  dependent  upon  the  optical  properties 
along  the  cloud  free  paths  of  sight  below  the  overcast. 

A  comparison  of  target  detection  distance  as  calcu¬ 
lated  with  and  without  an  overcast  cloud  are  shown  in  Fig. 
22.  The  diagnostic  calculations  again  were  made  using 
both  the  70-percent  and  99-percent  humidity  LOWTRAN 
rural  models  for  the  boundary  layer.  For  the  cloud  free 
cases,  calculations  are  shown  for  solar  zenith  angles  of  25, 
65  and  85  degrees.  A  solar  zenith  angle  of  55  degrees  and 
stratocumulus  overcast  of  average  optical  depth  were 
assumed  for  the  calculation  of  the  cloud  case.  As  dis¬ 
cussed  above,  the  slant  path  contrast  transmittance  and 
detection  distance  for  the  overcast  condition  are  not  sensi¬ 
tive  to  cloud  optical  depth  (or  solar  zenith  angle).  Notice 
in  Fig.  22  that  the  calculated  detection  distance  for  the 
cloud  case  corresponds  to  the  case  where  the  cloud  is 
removed  and  the  solar  zenith  angle  is  near  75  degrees. 
This  is  a  typical  result  for  conditions  where  a  Lambertian 
surface  reflectance  is  assumed  and  the  optical  properties 
along  the  viewing  path  are  the  same  for  both  the  clear  and 
cloud  examples. 

7.0  SUMMARY  COMMENTS 

The  analytic  techniques  that  were  developed  as  part 
of  the  experimental  optical  measurement  and  analysis  pro¬ 


gram  provide  relatively  fast  and  consistent  estimates  of  the 
visible  epectrum  contrast  transmittance  along  vertically 
inclined  paths  of  sight.  Essential  elements  of  the  scheme 
include  modeling  techniques  suitable  for  field  use  for  Ihe 
determination  of  the  apparent  background  sky,  cloud  and 
terrain  radiance  distribution  as  well  as  the  radiance 
transmittance  along  the  viewing  path.  Emphasis  was 
placed  upon  the  development  of  algorithms  which  are  sim¬ 
ple  and  easy  to  apply  yet  are  adaptable  so  as  to  lake  advan¬ 
tage  of  all  relevant  observations  and  forecast  information. 
The  procedures  developed  for  the  calculation  of  clear  sky 
radiance  distribution  are  directly  applicable  to  the  overcast 
sky  condition  through  introduction  of  the  appropriate  opti¬ 
cal  properties  of  the  cloud  layer.  Comparisons  of  model 
calculations  with  experimental  aircraft  measurements  of 
cloud  radiance  indicate  that  use  of  average  optical  depths 
for  various  cloud  types  as  derived  from  the  SOLMET  data 
base  yields  a  good  first  approximation  to  the  actual  over¬ 
cast  cloud  radiance  distribution. 

Another  recent  extension  to  the  image  transmission 
modeling  procedures  is  the  addition  of  techniques  to  deal 
with  the  combined  influences  of  atmospheric  and  t  -get 
factors  in  order  to  determine  target  detection  range  for 
objects  located  at  the  surface  and  viewed  from  aloft. 
Although  the  results  described  in  this  report  pertain  to 
simple  targets  with  uniform  backgrounds,  the  diagnostic 
model  calculations  provide  insight  into  the  relative  impor¬ 
tance  of  the  uncertainties  involved  in  the  estimates  of  tar¬ 
get  detection  distance.  Not  surprisingly,  the  determina¬ 
tions  are  most  sensitive  to  uncertainties  in  the  optical 
depth  of  the  atmosphere  which  in  turn  impact  directly  all 
three  components  of  contrast  transmittance,  i.e.  the 
inherent  and  apparent  background  radiance  and  the  path 
transmittance.  The  sensitivity  calculations  also  confirm 
the  strong  dependence  of  target  detection  distance  on 
knowledge  of  target  and  background  reflectivity.  The 
model  capability  for  joint  assessment  of  target  and 
environmental  effects  on  target  detection  range  also 
reveals  the  interactive  effects  of  such  parameters  as  solar 
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zenith  angle,  surface  reflectivity,  target  size,  and  visual 
fixation  time  for  overcast  and  clear  sky  conditions. 

Efforts  are  underway  to  improve  the  methodology 
for  estimating  image  transmission  characteristics  by  relaxa¬ 
tion  of  some  of  the  constraints  which  limit  the  applicability 
and  accuracy  of  the  techniques.  It  is  important  in  particu¬ 
lar,  to  develop  techniques  to  exploit  expanding  knowledge 
of  the  directional  dependence  of  the  reflectivity  of  various 
types  of  ground  cover.  Work  is  also  underway  to  incor¬ 
porate  analytic  techniques  to  estimate  slant-path  contrast 
transmittance  and  target  detection  distance  for  inclined 
paths  of  sight  in  fractional  cloud  cover  conditions. 
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APPENDIX  A 


SUMMARY  OF  INPUT  DATA, 
COMPUTATION  STEPS  AND 
RESULTANT  DATA 


A.l  General  Input  Data 

A. 1.1  Number  of  atmospheric  layers 

A.  1.2  Representative  wavelength  (/zm) 

A.  1.3  Extraterrestrial  solar  (lunar)  irradiance 
( W/m2nm ) 

A.  1.4  Rayleigh  optical  depth  of  top  layer 

A.  1.5  Solar  zenith  angles  (deg.) 

A.  1.6  Zenith  viewing  angles  (deg.) 

A. 1.7  Average  reflectance  of  underlying  surface  and 

reflectance  of  background  in  the  vicinity  of  target 
if  different  than  the  average  reflectance. 

A.2  Input  Data  Each  Layer 

A. 2.1  Base  altitude  (km) 

A. 2. 2  Scattering  mixing  ratio 

A.2.3  Single  scattering  albedo 

A. 2.4  Single  scattering  phase  function  (normalized  to 

integral  of  one) 

A.3  Computation  Steps 

A.3.1  Compute  Rayleigh  atmosphere  optical  thickness 
for  each  layer 

A.3. 2  Compute  total  optical  thickness  for  each  layer. 

A.3.3  Compute  optical  depth  for  the  base  of  each  layer. 

A.3.4  Compute  single  scattering  phase  function  and 
asymmetry  factors  for  each  layer.  (These  param¬ 
eters  may  be  specified  if  desired  in  lieu  of  the  use 
of  model  algorithms). 

A. 3. 5  Transform  optical  depth,  single  scattering  albedo, 

and  asymmetry  factor  in  accordance  with  delta- 
Eddington  approximation. 

A.3.6  Compute  path  function  distribution  for  the  base 
of  each  sub-layer  (each  layer  is  divided  equally 
into  2  sub-layers). 

A.3.7  Compute  path  radiance  distribution  for  each  sub¬ 
layer. 

A.3.8  Compute  path  transmittance  distribution  for  each 
sub-layer. 

A.3.9  Compute  sky/terrain  radiance  distribution  at 
sub-layer  boundaries. 


A.3. 10  Compute  slant  path  contrast  transmittance  of 
selected  paths. 


A. 4  Output  Tables  and  Data 

A. 4.1  Path  radiance,  L,  (z,9fi),  for  O°^0<85°  (upward 

looking  paths)  for  selected  levels  (max  5)  for 
each  0,  (i.e.  Sky  Radiance)  (IP/m^m  sr) 

A. 4. 2  Path  radiance.  Lr(z,9J3)  for  95°<0$18O° 
(downward  looking  paths)  for  selected  levels 
(max  5)  for  each  0,  ( W/m^m  sr) 

A. 4. 3  Path  plus  terrain  radiance  L,  (z,0j3)  for 

9500<  180°  for  selected  levels  (max  5)  for  each 
9,  sr) 

A. 4. 4  Contrast  transmittance,  Tc(zfi£)  between  object 

and  sensor  altitudes  (max.  5  altitude  intervals) 
for  0°O<  85°  and  for  each  0, . 

A. 4. 5  Contrast  transmittance  Tc(z,9Ji)  between  object 

and  sensor  altitudes  (max.  5)  for  95 °$0$  180° 
and  for  each  0, . 


A. 5  Additional  Output  From  The  Supplemental 
Eddington  Computer  Program 

The  standard  Eddington  computer  program  is 
employed  as  an  integrated  part  of  the  model  calculations 
of  contrast  transmittance.  The  following  output  from  the 
program  as  modified  by  the  delta-Eddington  transforma¬ 
tion  of  input  parameters  is  available  in  tabular  form  for 
the  base  of  each  sub-layer  and  for  each  solar  zenith  angle: 

•  optical  depths,  r  and  r' 

•  components  of  Eddington  diffuse  radiance,  LD  (: )  and 
Ld  (z )  ( W/m2fim  sr) 

•  total  scalar  irradiance  (H //m2fxm) 

•  downward  diffuse  irradiance  (W/mlnm) 

•  upward  diffuse  irradiance  (W!m1nm) 

•  solar  irradiance  (W/m2pm) 

•  total  downward  irradiance  ()V/m2fim) 
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APPENDIX  C 


GLOSSARY  AND  NOTATION 


The  notation  used  in  reports  and  journal  articles  produced  by  the  Visibility  Laboratory  staff  follows,  in  general,  the 
rules  set  forth  in  pages  499  and  S00,  Duntley  et  al.  (19S7).  These  rules  are: 

•  Each  optical  property  is  indicated  by  a  basic  (parent)  symbol. 

•  A  presubscript  may  be  used  with  the  parent  symbol  as  an  identifier,  e.g.,  b  indicates  background  while  t  denotes 
an  object,  i.e.  target. 

•  A  postsubscript  may  be  used  to  indicate  the  length  of  a  path  of  sight,  e.g.,  r  denotes  an  apparent  property  as 
measured  at  the  end  of  a  path  of  sight  of  length  r,  while  o  denotes  an  inherent  property  based  on  the  hypothetical 
concept  of  a  photometer  located  at  zero  distance  from  an  object,  i.e.  target. 

•  A  postsuperscript*  or  postsubscript,,  is  employed  as  a  mnemonic  symbol  signifying  that  the  radiometric  quantity 
has  been  generated  within  the  path  or  path  segment  by  the  scattering  of  ambient  light  reaching  the  path  from  all 
directions  and/or  by  emission. 

•  The  parenthetical  attachments  to  the  parent  symbol  denote  altitude  and  direction.  The  letter  z  indicates  altitude  in 
general;  z,  is  used  to  specify  the  altitude  of  a  target.  The  direction  of  a  path  of  sight  is  specified  by  the  zenith 
angle  0  and  the  azimuth  <t>  In  the  case  of  irradiances,  the  downwelling  irradiance  is  designated  by  d.  the 
upwelling  by  u. 

•  The  radiometric  symbols  used  herein  now  correspond  to  the  OSA  recommendations  in  Section  1  of  Driscoll  and 
Vaughn  (1978).  Prior  to  June  1980,  the  symbol  used  for  radiance  L  was  N,  for  irradiance  £  was  H,  and  for 
attenuation  length  £  was  L . 


Symbol  Units 


Quantity 


A(z) 

none 

Albedo  at  altitude  z,  defined 

A  (z)  -  £(z,u)  /  E(z,d). 

a(z ) 

m~' 

Absorption  coefficient. 

C(z) 

sr~l  m~' 

Diffuse  component  of  the  solar  almucantar  sky  radiance  ratio  to  the  sun  scalar 
irradiance  and  relative  airmass 

C(z)  —  J L(.z,9',<b')<r(.z,fi')d(l  /  se(z)  +  £»o(z)  /  st(z). 

c 

m/s 

Speed  of  light  c  —  2.99792458£8. 

C| 

Wm2 

First  radiation  constant  C|  -  In  he1  -  3.741832£-I6. 

Cj 

mK 

Second  radiation  constant  cj  —  chik  —  1.438786E-2. 

DM 

none 

Radiance  distribution  function  D(z)  -  J L(z,0',d>')sec9'dn  /  *(z) 

DM 

none 

Limb  darkening  factor  relating  the  average  sun  radiance  and  the  center  sun 
radiance  Z>(A)  -  tL0  /  SL0. 

W/m2 

Spectral  irradiance  (formerly  symbol  H)  defined  as  £t  —  j £„(z,0,d)  cos#  d fl 
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Symbol  Units  Quantity 


E 

W/m2  jxm 

oo 

Broad  band  sensor  irradiance,  defined  as  E  -  J" £»  Sx  7\  dk  /  8 k. 

0 

E(z,d) 

W/m 2 

Irradiance  produced  by  downwelling  flux  as  determined  on  a  horizontal  flat  plate 
at  altitude  z  (formerly  H(z,d)].  In  this  report  d  is  used  in  place  of  the  minus 
sign  in  the  notation  E//(z, — )]  which  appears  in  Duntley  (1969).  This  property 
may  be  defined  by  the  equation  E(z,d)  -  j L(z,9\<t>')  cos®' dd. 

E(z,u) 

W/m2 

Irradiance  produced  by  the  upwelling  flux  as  determined  on  a  horizontal  flat  plane 
at  altitude  z  [formerly  H(z,u)].  Here  u  is  substituted  for  the  plus  sign  formerly 
used  in  the  notation  l//(z,+)]. 

Hiz) 

m 

Scale  height  at  altitude  z.  the  height  of  a  homogeneous  atmosphere  having  the 
density  of  the  layer  at  altitude  z. 

h 

J/s 

Planck  constant  h  =  6.6261 76£" — 34. 

k 

J/K 

Boltzmann  constant  k  =  1 .38066267—23. 

lk 

W/sr  m2 

Spectral  radiance  (former  symbol  N). 

L 

W/sr  m2nm 

Broad  band  sensor  radiance  is  defined  as  L  =  J LKSKTK  dk  /  8\. 

0 

L0(z,.9,<b ) 

W/sr  m1 

Inherent  radiance  based  on  the  hypothetical  concept  of  a  photometer  located  at 
zero  distance  from  an  object  at  altitude  :  in  the  direction  specified  by  zenith  angle 

9  and  azimuth  <t>. 

L,(z,  9,<t>) 

W/sr  m2 

' 

Apparent  radiance  as  determined  at  altitude  z,  from  the  end  of  a  path  of  sight  of 
length  r  at  zenith  angle  9  and  azimuth  d>.  This  property  may  be  defined  by 

Lr(z,9,d> )  =  L0(z,.9,<t>)  Tr(z,9)  +  L’(z.9,d>). 

L  (k,  T) 

W/sr  m2 

Black  body  radiance  at  wavelength  k  and  temperature  T. 

; 

Lqiz,9,<b) 

W/sr  m 2 

Equilibrium  radiance  at  altitude  z  with  the  direction  of  the  path  of  sight  specified 

by  zenith  angle  9  and  azimuth  <t>.  This  property  is  a  point  function  of  position  ' 

and  direction. 

k 

j 

As  discussed  by  Duntley  ei  al.  (1957),  many  image  transmission  phenomena  are  i 

most  clearly  understood  in  terms  of  the  concept  of  equilibrium  radiance.  This 
concept  is  a  natural  consequence  of  the  equation  of  transfer  which  states  • 

analytically  that  in  any  path  segment  the  difference  between  the  output  and  input  j 

radiances  is  attributable  to  a  gain  term  and  a  loss  term,  such  that  some  unique  1 

equilibrium  radiance  Lq(z,9,<t>)  must  exist  at  each  point  such  that  the  loss  of  1 

radiance  within  the  path  segment  is  exactly  balanced  by  the  gain,  i.e.,  I 

M(z,9,d>)-  0.  I 

L 
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Symbol  Units 


Quantity 


By  virtue  of  this  concept  and  the  equation  of  transfer,  one  can  show  that  each 
segment  of  every  path  of  sight  has  associated  with  it  an  equilibrium  radiance,  and 
that  the  space  rate  of  change  in  image  forming  radiance  caused  by  the  path 
segment  is  in  such  a  direction  as  to  cause  the  output  radiance  to  be  closer  to  the 
equilibrium  radiance  than  is  the  input  radiance.  This  segment  by  segment 
convergence  of  the  apparent  radiance  of  the  object  field  to  the  dynamic 
equilibrium  radiance  was  clearly  illustrated  by  the  data  in  the  1957  paper 
referenced  above. 


W/sr  m 2 

Path  function  at  altitude  z  with  the  direction  of  the  path  of  sight  specified  by 
zenith  angle  8  and  azimuth  </>.  This  property  is  defined  by  the  equation 
L,(z,8,</>)  =  L.s(z,8,<t>)  +  Z.,fl(z). 

L.s(z,8,<t>) 

W/sr  m2 

Path  function  due  to  scattering,  defined  by  L.(z,0,<£)  =  f  cr(z,0')L(z.0',<i>')d(l' 

An 

L.a(z) 

W/sr  m 2 

Path  function  due  to  emittance  (absorption)  L,a(z )  =  a(z)  L(k,T). 

L,’(z,e,0) 

W/sr  m2 

Path  radiance  as  determined  at  altitude  z  at  the  end  of  a  path  of  sight  of  length  r 
in  the  direction  of  zenith  angle  8  and  azimuth  <£. 

Ll(z,8,<f>) 

W/sr  m2 

Sky  radiance  at  altitude  z,  zenith  angle  8  and  azimuth  <t>  Also  the  path  radiance 
for  the  path  of  sight  of  length  °°  from  out  of  the  atmosphere  to  altitude  z. 

,L„(z,0 

W/sr  m2 

Apparent  radiance  of  the  center  of  the  solar  disk  as  determined  from  the  end  of  a 
path  of  sight  of  length  «=  from  out  of  the  atmosphere  to  altitude  z  at  the  zenith 
angle  of  the  sun 

£(z) 

m 

Attenuation  length  at  altitude  z.  The  attenuation  length  is  the  distance  at  which 
the  signal  is  attenuated  to  1/e. 

Af(X,T) 

W/m2 

Black  body  exitance  (emittance)  M(\,T)  -  nLik.T). 

W/m 2 

Black  body  scalar  exitance  (emittance)  m(k,T)  -  4trL(\<T). 

m.(z) 

W/m 3 

Scalar  exitance  (emittance)  per  length  m,(z)  -  / L.(z,8,<t>)d ft. 

4» 

m,s(z ) 

W/m 3 

Scalar  exitance  (emittance)  per  length  due  to  scattering 

m.5(z)  -  jTz..j(z,0,<£)</ft  -  «(z)s(z). 

W/m2 

Scalar  exitance  (emittance)  per  length  due  to  absorption  m,a(z)  -  a(z)mik.T). 

mm(.z,9 ) 

kf/m1 

o» 

Absolute  air  mass  at  angle  9  mm(z,9)  -  f  p(z)dr 

X 
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Symbol  Units  Quantity 

«o.(z,0)/m<,(z, 0°)  none  Relative  optical  air  mass. 

none  Relative  optica!  airmass  at  zenith  angle  of  the  sun,  shorthand  for 

mo.(z,@s)/m„,(z,0°). 

A'(z)  none  Refractive  modulus  l V(z)  *=  [n(z)  -  l]£6. 


n(z)  none  Refractive  index. 


@(z>  none  Optical  scattering  mixing  ratio  at  altitude  z.  This  quantity  is  defined  as  the  ratio 

of  the  total  volume  scattering  coefficient  at  altitude  z4  to  the  molecular  (or 
Rayleigh)  volume  scattering  coefficient  at  the  same  altitude  z. 
Q(z)  =  s(:)/ks(z). 

2  (z,0)  none  Volume  scattering  function  ratio  at  altitude  z.  This  quantity  is  defined  as  the  ratio 

of  the  total  volume  scattering  function  at  altitude  z  and  scattering  angle  0,  to  the 
molecular  (or  Rayleigh)  volume  scattering  function  at  the  same  altitude  and 
scattering  angle.  2(z,0)  -  o-(z,/ 3)  /  *o-(z,/ 3). 

r  m  Path  length,  for  paths  of  sight  at  zenith  angles  0  to  70  degrees,  r  =  secB  A z. 

T,  none  Standardized  relative  spectral  response  of  filter/cathode  combination  where  is 

spectral  sensitivity  of  the  multiplier  phototube  cathode  and  T*  is  spectral 
transmittance  of  optical  filter.  The  relative  spectral  response  values  are 
normalized  to  the  peak  value. 

stz)  m  1  Total  volume  scattering  coefficient  as  determined  at  altitude  z.  This  property  may 

be  defined  by  the  equations 

s(z)  =  J*c z(z,0)dft  =  ftS(z)  +  ms(z). 

Air 

In  the  absence  of  atmospheric  absorption,  the  total  volume  scattering  coefficient  is 
numerically  equal  to  the  attenuation  coefficient. 

Volume  scattering  coefficient  for  Mie  i.e.  particulate  or  droplet,  scattering  at 
altitude  z. 

rs(z)  m  '  Volume  scattering  coefficient  for  Rayleigh  i.e.  molecular,  scattering  at  altitude  z. 

T  K  Absolute  temperature  in  degrees  Kelvin. 


T,(z,9)  none  Radiance  transmittance  as  determined  at  altitude  z  for  a  path  of  sight  of  length  r 

at  zenith  angle  9  (formerly  referred  to  as  "beam"  transmittance).  This  property  is 

independent  of  azimuth  in  atmospheres  having  horizontal  uniformity.  It  is  always 
the  same  for  the  designated  path  of  sight  or  its  reciprocal. 

T»(r,9)  none  Tladiance  transmitance  for  the  path  of  sight  at  zenith  angle  9  from  out  of  the 

atmosphere  to  the  altitude  z. 
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Symbol 

Units 

Quantity 

none 

Optical  thickness  /i2(z)  -  J*  a(z)dz 

none 

Absorption  optical  thickness  atlz(z)  -  f  o(z)<Jlz. 

none 

Mie  optical  thickness  Mr42(z)  -  J"  Ms(z)dz. 

2, 

*<12^ 

none 

2 

Rayleigh  or  molecular  optical  thickness  Rt^(z)  -  J"  Rs(z)et. 

w 

J/m} 

Radiant  density. 

W 

W 

Watt  W  -  S/s. 

Z 

m 

Altitude,  usually  used  as  above  ground  level. 

h 

m 

Altitude  of  any  applicable  object  or  target 

a(z) 

m' 

Volume  attenuation  coefficient  as  determined  at  altitude  a(z)  =  a(z)  +  s<z). 

P 

deg 

Symbol  for  scattering  angle  of  flux  from  a  light  source.  It  is  equal  to  the  angle 

between  the  line  from  the  source  to  any  unit  scattering  volume  and  the  path  of  a 
ray  scattered  off  this  direct  line.  See  illustration. 


V 


A  none  Symbol  to  indicate  incremental  quantity  and  used  with  r  and  c  to  indicate  small. 

discrete  increments  in  path  length  r  and  altitude  : 

#X  nm  Effective  passband  (formerly  designated  "response  area"!  for  a  filtered  sensor  is 

defined  as  8A  —  J*  Sk  Tk  d\ 

0 
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Symbol 

Units 

Quantity 

e(z) 

W/m1 

Scalar  irradiance.  This  may  be  defined  as  the  radiant  flux  arriving  at  a  point,  from 
all  directions  about  that  point,  at  altitude  z  (Tyler  and  Preisendorfer  (1962)). 
e (z)  L(z,0',<f>')da'. 

d«(z) 

W/m1 

Diffuse  scalar  irradiance  rf«(z)  -  t(z)  -  se(z). 

5«(r) 

W/m2 

Sun  scalar  irradiance  at  altitude  z  se(z)  -  s*(°o)  T„(z,9s). 

( 

m 

Radius  of  the  earth. 

0 

deg 

Symbol  for  zenith  angle.  This  symbol  is  usually  used  as  one  of  two  coordinates  to 
specify  the  direction  of  a  path  of  sight. 

0‘ 

deg 

Symbol  for  zenith  angle  usually  used  as  one  of  two  coordinates  to  specify  the 
direction  of  a  discrete  portion  of  the  sky. 

»S 

deg 

Zenith  angle  of  the  sun. 

A 

nm 

Symbol  for  wavelength. 

A 

nm 

Mean  wavelength  is  defined  as  A  -  f  A  SKTk  dK  /  SA. 

0 

/x(z,/3) 

sr1 

Ratio  of  solar  almucantar  sky  radiance  to  scalar  sun  irradiance  and  airmass. 
p(z,/3)  -  Li,  (z,9s,fi)  /  (  s€(z)m(0j)]. 

W/mi 

Net  irradiance.  f(z)  -  E(z,d)  -  E(z,u). 

p(z) 

kg/  m} 

Density  at  altitude  z. 

tr(z,/3) 

m-'sr"1 

Symbol  for  volume  scattering  function.  Parenthetical  symbols  are  z  to  designate 
altitude  and  0  to  designate  the  scattering  angle  from  a  source. 

tr(z,/9)/s(z) 

sr~l 

Normalized  volume  scattering  function.  This  may  be  defined  by  the  equation 
£[  a(z,l 3)  /  s(z)]rffl  -  1. 

t4,(z.^) 

sr 

Optical  thickness  function.  r4z(z,/ 3)  -  J* cr(z,/3i<fc. 
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Symbol  Units 


Quantity 


* 

deg 

Symbol  for  azimuth.  The  azimuth  is  the  angle  in  the  horizontal  plane  of  the 
observer  between  a  fixed  point  and  the  path  of  sight.  The  fixed  point  may  be  for 
example,  true  North,  the  bearing  of  the  sun,  or  the  bearing  of  the  moon  This 
symbol  is  usually  used  as  one  of  two  coordinates  to  specify  the  direction  of  a  path 
of  sight. 

deg 

This  symbol  for  azimuth  is  usually  used  as  one  of  two  coordinates  to  specify  the 
direction  of  a  discrete  portion  of  the  sky. 

'll 

deg 

Angular  solar  radius  at  true  earth-to-sun  distance. 

* 

deg 

Angular  solar  radius  at  mean  solar  distance. 

n 

sr 

Symbol  for  solid  angle  For  a  hemisphere,  fl  -  2 n  steradians.  For  a  sphere, 
fl  -  4jt  steradians. 

<0(2) 

none 

Single  scattering  albedo  <u(z)  -  s(z)  /  a(z). 
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Table  C.l.  Notaiional  equivalencies 


Chandrasekhar  (1960) 

Visibility  Laboratory 

Symbol 

Definition 

Symbol 

Definition 

E 

Radiant  energy 

E 

Irradiance 

TT  F 

Net  flux 

f 

Net  irradiance 

1 

Specific  intensity 

H:.9.<b) 

Radiance 

J 

Average  intensity 

t/4rr 

Average  radiance 

j 

Mass  emission  coefficient 

Path  function/density 

p(z) 

Jp 

Mass  emission  coefficient  x  density 

L.(:,9,  <t>) 

Path  function 

,ls> 

Mass  emission  coefficient  due 

Path  function  due  to  scattering 

J 

to  scattering 

p(:) 

divided  by  the  density 

k 

Mass  absorption  coefficient 

alp 

Attenuation  coefficient  divided 
by  the  density 

kp 

Mass  absorption  coefficient  x  density 

a(z) 

Attenuation  coefficient 

p  ( ««(-)) 

Phase  function 

4n(r(p)/a 

Four  tt  times  volume  scattering 
function/attenuation  coefficient 

s 

Thickness 

r 

Path  length 

3 

Source  function 

Equilibrium  radiance 

U 

Integrated  energy  density 

w 

Radiant  density 

Single  scattering  albedo 

<u(z) 

Single  scattering  albedo 

H 

Scattering  angle 

p 

Scattering  angle 

9 

Polar  angle 

9 

Zenith  angle 

M 

cosine  of  polar  angle 

cos  9 

cosine  of  zenith  angle 

V 

Frequency 

1/A 

Inverse  wavelength 

p 

Density 

P 

Density 

T 

Normal  optical  thickness 

Optical  thickness 

V 

Azimuth 

<b 

Azimuth 

(V 

Solid  angle 

n 

Solid  angle 
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